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Praticamente todas as máquinas (2) 
equipamentos estão equipados com 
o-rings. Na construção mecânica, na 
industria electrotécnica, em mecânica 


de precisão, na indústria sanitária, na 


construção de acessórios e em muitos 
outros campos, os o-rings são ele- 
mentos de vedação indispensáveis. 
Os o-rings SIMRIT têm grandes van- 
tagens. 

O A qualidade. Ela resulta da eficiência 
dos nossos laboratórios, do interesse 
que dedicamos à investigação e do 
acabamento do material. No nosso 
trabalho são utilizados os mais per- 
feitos e recentes meios da mais mo- 
derna e avançada técnica. 

6 Grande variedade de medidas. Pro- 
duzimos o-rings em mais de 2.500 di- 
mensões e temos, permanentemente, 


O Grand 


de SIMRIT à basesde NBR, 
rezas 70 e 90 Shore A, e de Viton. 


E para casos especiais, utilizamos 


outros materiais. 
ham de suportar temper 
mas ou pressões elevadas, q 
destinem a vedar os produtos m 
agressivos, a verdade é que estamos 
sempre aktos a escolher o material 
exacto para o seu caso. 
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jectista. 
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movimento relativo — vedação dinã- 
mica. 
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Projecto apoiado em computador de sistemas de controlo 
Bloco operacional 


RESUMO 


O artigo descreve com algum detalhe o bloco ope- 
racional de um conjunto articulado de programas para o 
projecto apoiado em computador de sistemas de controlo, 
Neste bloco encontram-se implementados os métodos do 
lugar geométrico das raizes, Nyquist, Bode e da simula- 
ção digital da resposta no tempo. Apresenta-se os algo- 
ritmos utilizados bem como alguns exemplos de aplicação. 


0 — INTRODUÇÃO 


Antes de passarmos à análise do bloco operacional 
de um cenjunto articulado de programas para o projecto 
apoiado em computador de sistemas de controlo, convém 
referir que este artigo se insere e é continuação de um 
outro artigo já publicado [1] onde se descreve o bloco 
de manipulação e tratamento de dados, sendo ainda de 
frisar que o presente trabalho resulta de uma colaboração 
entre o Laboratório Nacional de Engenharia Civil (L, N. 
E. C) e o Centro de Automática da Universidade Técnica 
du Lisboa (C, A. U, T, L). 

No bloco operacional do referido conjunto articulado 
de programas, <CONPAC>S, foram implementadas as 
técnicas ou métodos mais vulgarmente utilizados na aná- 
lise e sintese de sistemas de controlo univariáveis (uma 
entrada/uma saida) e que são os seguintes: 


— Método do lugar geométrico das raízes («root 
locus»), 

— Método de Nyquist, 

— Método de Bode, 


— Simulação digital da resposta no tempo, 


sendo de referir que em todos eles se admitem sistemas 
que apresentem atraso no tempo. 


Como é sabido, os resultados obtidos com estas 
técnicas exprimem-se, fundamentalmente, através de dia- 


(7 Especialista do Laboratório Nacional de Engenharia Civil. 


FRANCISCO CARVALHAL (') 
ISABEL LOURTIE (7) 
PEDRO LOURTIE () 


ABSTRACT 


The paper describes with some detail the operational 
block of a control systems computer-aided design packa- 
ge. The root locus, Nyquist, Bode and time response 
methods are implemented in this block. The algorithms 
are presented as well as some application examples. 


gramas; sendo assim, a parte respeitante à apresentação 
gráfica dos resultados assume especial relevância, Este 
facto foi, com efeito, tomado em devida consideração na 
elaboração e estruturação dos algoritmos adoptados. Na 
realidade, como adiante se verá, os pontos dos diagramas 
obtidos para qualquer uma das técnicas referidas estão 
graficamente equidistantes entre si, o que melhora consi- 
deravelmente a observação daqueles diagramas. 

Na exposição dos algoritmos utilizados, tentaremos, 
essencialmente, descrever a «filosofia» que preside ao 
seu funcionamento, evitando-se sobrecarregá-la com ex- 
cessivos detalhes da implementação, detalhes estes que 
poderiam, eventualmente, desmotivar o leitor, Reportar- 
-nos-emos a sistemas cujo diagrama de blocos simplifi- 
cado se representa na Fig. 1.1. 


FIG. 1.1 
G — sistema controlado 


C — controlador 
H — retrúacção 


() Assistente do Instituto Superior Técnico, Investigador do Centro de Automática da Universidade Técnica de Lisboa (]. N.1. C.). 
() Professor do Instituto Superior Técnico, Investigador do Centro de Automática da Universidade Técnica de Lisboa (Il. N. 1. C.). 


Manuscrito recebido para publicação em 3/7/50. 
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0.1 — Entrada no bloco operacional 


A entrada no bloco operacional para a análise de um 
sistema pressupõe que já houve, anteriormente, uma ou 
mais passagens pelo bloco de manipulação de dados de 
modo a ter-se possibilitado a criação dos ficheiros neces- 
sários, contendo a informação requerida, adequadamente 
organizada [1]. 

Nestas circunstâncias, qualquer que seja a técnica 
seleccionada, é obrigatório seguir-se um procedimento 
preparatório, o qual envolve algumas perguntas dirigidas 
ao utilizador. Assim, é indagado, em primeiro lugar, se se 
pretende a análise do sistema em anel fechado ou em 
anel aberto, No primeiro caso, é lido o ficheiro que con- 
tém a função de transferência do sistema em anel fecha- 
do entrando-se, em seguida, no método escolhido, No 
segundo caso, sistema em anel aberto, é lido o ficheiro 
que contém a informação relativa a G(s) e, de seguida, 
é averiguado se já foram criados os ficheiros corres- 
pondentes aos sub-sistemas controlador C(s) e de re- 
troacção H(s) ou, alternativamente, se se pretende consi- 
derar C(s)=1 e/ou H(s) = 1; finalmente, o produto 
G(s).C(s).H(s) é calculado antes de se entrar na técnica 
seleccionada. Em qualquer destes dois casos, anel fecha- 
do ou anel aberto, o procedimento referido provoca um 
retorno automático ao bloco supervisor principal, prece- 
dido de uma mensagem apropriada para o utilizador, se 
não tiverem sido já criados todos os ficheiros necessá- 
rios. 


1 — LUGAR GEOMÉTRICO DAS RAÍZES 
(«ROOT LOCUS-») 


O algoritmo elaborado para o traçado do diagrama 
do lugar geométrico das razes (LGR) recorre, fundamen- 
talmente, a uma busca sistemática, sequencial e angular, 
dos pontos que constituem os diferentes ramos do LGR 
de um sistema, para ganhos positivos, com excepção 
dos ramos (ou de troços a eles pertencentes) que coin- 
cidam com o eixo real (abcissas), os quais são traçados 
utilizando a propriedade referida em 1.2.1 [2, 3]. 

A entrada do programa relativo a esta técnica, é 
conhecida a função de transferência do sistema a analisar 
quer na sua forma racional quer na sua forma factorizada, 
ou seja, são conhecidos os valores não só dos coeficien- 
tes dos polinómios do numerador e do denominador como 
também dos polos e dos zeros, valores estes armazena- 
dos em variáveis indexadas ou «arrays» adequados; é 
conhecido ainda o valor de «x, respeitante ao atraso no 
tempo, o qual pode ser nulo (sistema sem atraso). Antes 
do traçado do LGR propriamente dito, são executadas 
diversas operações preliminares; durante esta fase inicial 
o programa requer ainda alguns dados fornecidos pelo 
utilizador, como adiante se verá. Na fase final, o progra- 
ma permite ainda, caso se pretenda, não só o traçado 
de linhas (rectas) de amortecimento relativo constante 
como também a localização dos polos do sistema (com 
o anel de retroacção fechado) para diversos valores do 
ganho em anel aberto, 


1.1 — Fase preliminar ou inicial 


Na fase inicial já atrás referida, o programa executa 
as acções e operações que a seguir se indicam: 


1) Dá o valor zero a determinadas variáveis inteiras 
com funções do tipo «flag». 

2) Pergunta ao utilizador quais as dimensões exte- 
riores (em cm) do desenho correspondente ao 
LGR. 

3) Abre o ficheiro (de extensão PLT) que conterá 
toda a informação que permitirá, posteriormente, 
representar-se o diagrama quer no «plotter» quer 
num terminal gráfico. Nota: este procedimento 
adapta-se à configuração actual do computador 
digital DEC-10 do LNEC. 

4) Coloca em três «arrays» diferentes os polos 
complexos, os zeros complexos e os polos e 
zeros reais. Em cada um deles, os polos e os 
zeros ficam ordenados por ordem decrescente 
dos valores das respectivas partes reais. Aten- 
dendo-se à simetria existente relativamente ao 
eixo real, não se incluem neste arranjo os polos 
e zeros conjugados de parte imaginária negativa. 
Com efeito, devido à referida simetria, todas as 
buscas são realizadas apenas no semi-plano cor- 
respondente à parte imaginária positiva, 

5) Calcula os valores extremos de ambos os eixos 
de modo a definir uma área minima que inclua 
todos os polos e zeros. Indica estes valores ao 
utilizador permitindo que este os modifique ou 
os mantenha. 

6) Calcula o valor para o incremento, ou passo, a 
utilizar na busca dos pontos dos ramos do LGR. 
O critério utilizado baseia-se na dimensão da 
área definida e nas distâncias dos polos comple- 
xos aos zeros complexos e fronteira daquela 
área. Indica o valor calculado ao utilizador, per- 
mitindo também que este o modifique caso o 
pretenda. 

7) Traça e gradua os eixos. Representa todos os 
polos por uma cruz (x) e todos os zeros por um 
losango (9). 

8) Determina quais os pontos singulares, 


1.1.1 — Cálculo dos pontos singulares 


A função de transferência racional de um sistema 
em anel aberto pode ser representada como segue 


K.N(s) =as 
FE) = mas 6 (1.1) 


Ds) 
em que 
F(s) = C(s). G(s). H(s) 


e, por conseguinte, 


N(s) = Ne(s) . Na(s) . Ny(s) 
D(s) = De(s). Da(s) . Dy(s) 
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e em que, por outro lado, & representa o valor do atraso 
no tempo e K o valor do ganho; s = 7 + ju é a variável 
complexa. 

Um ponto s* qualquer que satisfaça a relação 


Esp (1.2) 


é designado por ponto singular [4]. Atendendo à expres- 
são (1.1), da expressão (1.2) resulta que 


N(s). D'(s) — N'(s). D(s) + aN(s). D(s) = 0 (1.3) 


s a 
dado que el +*0 e que N(s) & 0 para valores finitos 
do ganho K. Se não se verificar atraso no tempo, & = 0, 
a expressão anterior devem 


N(s) . D'(s) — N'(s). D(s) = 0 (1.4) 


Designando por m a ordem do polinômio N(s) e por 
na ordem de D(s) verifica-se que existirão (n+m) pon- 
tos singulares para «-£0 e (n+m—1) pontos singulares 
sa 4=0. 

No programa a determinação dos pontos singulares 
consiste, como é natural, no cálculo das raizes dos poli- 
nómios representados em (1.3) ou (1.4), o que se faz 
recorrendo, sucessivamente às subrotinas PDER (deri- 
vação de polinómios), PMPY (multiplicação de polinó- 
mios) e PRBM (determinação de raizes reais e complexas 
de polinómios, com utilização da iteração de Bairstow). 
Estas subrotinas pertencem à Scientific Subroutine Packa- 
ge (55P) da IBM, residente no computador DEC-IO do 
LNEC. 


Para cada um dos pontos singulares assim obtidos 
é calculado o valor correspondente do ganho K, sendo 
eliminados aqueles para os quais se constata que K < 0. 

A importância da determinação dos pontos singula- 
res consiste não só no facto de neles se incluirem os 
pontos de saida e de entrada de ramos no eixo real 
(pontos «break-away» e «break-in»), pontos singulares 
reais, mas também por indicarem pontos onde dois ou 
mais ramos distintos do LGAR se intersectam ou tocam. 


1.1.2 — Valor do Ganho 


Para um sistema com retroacção negativa, a função 
de transferência em anel fechado é dada por 


K. NASINA(S)D ji 
T(s) = Pp ca 4 (1.5) 
D(s)+K.N(s)e - 


atendendo a (1.1) e expressões subsequentes, Ássim, 
para diversos valores do ganho K, os polos do sistema 


em anel fechado são dados pelas raizes ou zeros da 
função complexa 


a(s) = D(s) + K.N(s) e “ =0 (1.6) 
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ou, caso não exista atraso no tempo, 
g(s) = D(s) + K.N(s) = 0 (1.7) 


Tem-se, portanto, que o ganho K em cada ponto s 
do LGR obedece à condição expressa por 


K = E (1.8) 


N(s) e É 


Tomando a parte real e a parte maginária da expres- 
são (1.6) verifica-se que o valor do ganho, para um 
ponto s qualquer pertencente ao LGR, é dado [4] por 


Re(D(s)) 
K = — MDA (1.9) 


Re(N(s) e a 
ou, 


Im(D(s)) 


cs (1.10) 
Im(N(s) e dá , 


FÉ ras 


1.2 — Cálculo e traçado dos ramos 


Concluída a fase preliminar atrás descrita, inicia-se 
a execução do traçado propriamente dito do LGR, ou 
seja, o cálculo e traçado dos diferentes ramos; como é 
sabido, o número de ramos distintos é igual ao número 
de polos da função de transferência em anel aberto, 
quando não existe atraso no tempo. O procedimento 
sequencial utilizado é o que a seguir se indica: 


1) Traça os ramos existentes no eixo real, os quais 
designaremos por ramos reais. 

2) Traça os ramos que eventualmente intersectem a 
fronteira da área de interesse, previamente defi- 
nida pelo utilizador. Como foi referido na parte 
inicial de 1.1, esta área (e, portanto, a fronteira 
que a limita) refere-se apenas ao semi-plano cor- 
respondente à parte imaginária positiva, 

3) Traça os ramos que partem dos polos complexos, 
Eventualmente, se um ramo traçado a partir da 
fronteira atingiu um destes polos, não é então 
aqui considerado uma vez que, como é evidente, 
o ramo respectivo já se encontra traçado, 


4) Se, por algum motivo, algum ou alguns dos zeros 
complexos, caso existam, não forem atingidos por 
nenhum ramo, inicia-se então o seu traçado a 
partir destes zeros complexos, 


Com excepção dos ramos reais, os pontos dos res- 
tantes ramos, que designaremos por ramos complexos, 
são determinados através de um mesmo procedimento, 
o qual se baseia, essencialmente, nos valores assumidos 
pelo que denominamos de função de desvio angular 
construída com base na condição de ângulo ou condição 
angular. 
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1.2.1 — Ramos reais 


O traçado dos ramos reais é feito com base na pro- 
priedade do LGR que refere que para sistemas com 
retroacção negativa e valores positivos do ganho em anel 
aberto estes ramos existem nos segmentos do eixo real 
para os quais a soma de polos e zeros à sua direita 
(maiores valores de x) é um número impar [2]. 


1.2.2 — Condição angular, Função de desvio angular 


Se um ponto qualquer s=x+iy do plano pertence 
a um ramo do LGR deve então satisfazer a condição 
angular, a qual se pode representar pela expressão se- 
quinte 


[r= 001,0, 1,..] (141) 


em que Y; e O, são ângulos ou contribuições angulares 
referentes aos zeros e polos, respectivamente. O signifi- 
cado destas contribuições angulares pode ser facilmente 
apreendido pela observação da Fig. 1.2, na qual se apre- 
senta um exemplo do cálculo de Q para um ponto s = 
= x+ijyv. 


FIG. 1.2 


A busca dos pontos dos ramos complexos baseia-se 
na referida condição angular ou, mais concretamente, 
numa função de £)? procurando detectar, como adiante 
se verá, dois valores daquela função com sinais contrá- 
rios e de tal modo que para os pontos que satisfaçam 
aquela condição a função assuma o valor zero, 

Observando-se a expressão (1.11) verifica-se que a 
condição angular é satisfeita quando £? toma o conjunto 
de valores (..., —540º, —180º, 180º, 540º ...), apresen- 
tando, por conseguinte, uma periodicidade igual a 360º. 

Tendo em atenção o que se expôs nos parágrafos 
anteriores, optou-se pela seguinte função, 


F(O) = cos (0/2) (1.12) 


que designamos, como já foi atrás referido, por função 
de desvio angular, 

Dado que esta função tem como variável indepen- 
dente um somatório dos ângulos Y, e O, os quais, por 
sua vez, são calculados a partir da função circular inversa 
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arc tan, tornam-se necessárias algumas precauções, as 
quais se encontram implementadas no programa, precau- 
ções estas que se podem considerar como normais quan- 
do se lida com ângulos e funções circulares inversas. 


1.2.3 — Pesquisa dos ramos complexos 


Cada ponto de um ramo complexo é determinado 
através de uma pesquisa circular, com base nos valores 
assumidos pela função desvio é tendo em conta os dois 
pontos anteriormente calculados, O e P, têm coordenadas 
(x-1, Y-1) € (x;, y;), respectivamente, A pesquisa do 
novo ponto faz-se a partir do ponto P, o último a ser 
determinado. Em primeiro lugar, calcula-se o ângulo, [Á, 
que o segmento OP faz com o eixo real, com base no 
conhecimento das coordenadas atrás referidas, para a 
determinação do primeiro ponto, feita quando se inicia 
o algoritmo de pesquisa do ramo, é dado um valor de [À 
que depende do ponto de partida, À pesquisa do novo 


5 
si 
| “og 18 
ôs Fá a 
Favo ; nE- 
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Ls 
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FIG. 1.3 


ponto do ramo efectua-se sobre os pontos de uma cir- 
cunferência cujo raio é As, valor do incremento ou passo 
já referido em 1.1, e com centro em P. Um primeiro 
ponto de pesquisa, S,, é então determinado no prolon- 
gamento do segmento OP, com coordenadas dadas, por- 
tanto, pelas expressões 


Xo =X; + As.cos ([Á) 
(1.13) 
Yo =y + As.sen (|) 


Para este ponto S, são então calculados, através 
das expressões (1.11) e (1.12), os valores de 9º e de 
H(0)). Caso este último valor seja nulo ou, mais corree- 
tamente, caso 


HO) | < Af (1.14) 


o ponto S, é considerado como o novo ponto; Af é um 
valor pequeno e constante, internamento fixado. Caso 
não se verifique (1.14) inicia-se então uma busca sequen- 
cial para pontos alternadamente à direita e à esquerda 
de 3, (pontos Sy € Sj.). Mais concretamente, a sequên- 
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cia de pontos é constituida por Sig, Sjer S2g, Sã. , etc. 
As coordenadas de cada um destes pontos de pesquisa 
são calculadas por expressões semelhantes às (1.13) em 
que os ângulos são dados por 


Pua=B-k.aú 
Nr = À +k.am 


sendo A/Á um incremento angular constante e interna- 
mente estabelecido. Para cada um destes pontos deter- 
mina-se o valor da função desvio, isto é, fl) ou 
fr(O). Logo que se detecta uma inversão de sinal da 
função desvio entre dois pontos consecutivos, entra-se 
numa rotina onde se determina o ponto pertencente ao 
ramo (para o qual a expressão (1.14) é satisfeita) por 
intermédio de um algoritmo baseado no método da falsa 
posição [5], o qual apresenta a caracteristica de ser 
sempre convergente. 


Após ter sido determinado o novo ponto, entra-se 
numa rotina de verificação que averigua se aquele se 
encontra dentro de uma vizinhança de raio As do eixo 
real, da fronteira, de um ponto singular complexo, de um 
zero complexo ou, finalmente, de um polo complexo. 
Caso se verifique a ocorrência de qualquer um destes 
casos, a pesquisa de novos pontos do ramo cessa; caso 
contrário, continua. Relativamente aos polos e zeros com- 
plexos, aquela ocorrência implica o levantamento de uma 
«flag» correspondente, indicando que aquele polo (ou 
zero) foi já atingido por um ramo, Deve referir-se que 
esta rotina de verificação não é aplicada aos dois pri- 
meiros pontos determinados, atendendo ao critério que 
define o valor do incremento As. 


Convém ainda ser referido que o algoritmo adoptado 
permite que os pontos dos ramos do LGR sejam equi- 
distantes entre si, Por outro lado, a experiência tem 
demonstrado que é de execução rápida, 


1.2.4 — Pesquisa da fronteira 


O algoritmo da pesquisa da fronteira da área de 
interesse é muito semelhante ao algoritmo referente à 
busca dos ramos complexos. Como já atrás se referiu, 
esta pesquisa destina-se a investigar a existência de 
ramos que, eventualmente, intersectem a fronteira; deter- 
minado um ponto de intersecção, recorre-se então ao 
algoritmo de traçado de ramos complexos, O objectivo 
primário desta pesquisa consiste em tornar possível o 
traçado de ramos que, em determinadas circunstâncias, 
não o seriam; incluem-se nestas circunstâncias, factores 
tais como dimensões da área de interesse e/ou existência 
de atraso no tempo, como adiante se verá num dos 
exemplos de aplicação apresentados. 


Para determinar as coordenadas dos eventuais pon- 
tos de intersecção, a fronteira é dividida em três zonas: 
duas em que y varia desde As até Ymax & X toma, res- 
pectivamente, os valores constantes Xin € Xmaw € UMa 
em que x varia desde Xin até Xyax COM Y constante e 
igual a Ymax- Conforme os valores constantes das coor- 
denadas dos pontos de pesquisa sejam Xin. Ymax OU Xmax 
os ângulos iniciais para a pesquisa dos ramos complexos 
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a partir dos pontos de intersecção são iguais a 0º, 270º 
e 180º, respectivamente. 

Em cada uma das três zonas da fronteira são consi- 
derados diversos pontos equidistantes entre si de As. 
Sequencialmente, verifica-se para cada um deles se já 
foi atingido previamente por um ramo (para tal existe um 
«array» adequado que regista as coordenadas de pontos 
da fronteira já atingidos por ramos). Caso tal se verifique, 
o ponto não é considerado; caso contrário, calcula-se o 
valor da função de desvio angular correspondente a esse 
ponto, Quando se detecta uma inversão do sinal daquela 
função entre dois pontos consecutivos, aplica-se então 
o método da falsa posição que dará as coordenadas do 
ponto de intersecção, para o qual à condição (1.14) de- 
verá ser satisfeita. A pesquisa da fronteira só termina 
quando todos os pontos sobre ela considerados tiverem 
sido investigados. 


1.2.5 — Ramos a partir de polos e zeros complexos 


Para cada polo e zero complexo verifica-se, pelo 
estado da «flag» respectiva já atrás referida, se um ramo 
ja os atingiu. Em caso afirmativo não são considerados. 
Caso contrário, determina-se o ângulo de partida (polo) 
ou ângulo de chegada (zero) [2], cujo valor é fornecido 
à rotina de traçado de ramos complexos, a qual é então 
activada, 


1.3 — Fase final 


Concluído o traçado de todos os ramos do LGR, o 
utilizador pode ainda optar pelo traçado de linhas de 
amortecimento constante (para valores desse amorteci- 
mento não nulos e inferior a 1) e pela localização dos 
polos do sistema em anel fechado para diversos valores 
do ganho em anel aberto, 

As linhas de amortecimento relativo constante são 
traçadas sabendo-se que se tratam de rectas, passando 
pela origem, e fazendo um ângulo com o semi-eixo nega- 
tivo de x dado por arc cos É, em que É representa o 
valor do amortecimento relativo [2]. 

A localização dos polos do sistema em anel fechado 
é realizada através da determinação dos zeros da função 
(1.6) ou, caso não exista atraso no tempo, de (1.7). 

Para a determinação dos zeros de (1.6) utiliza-se a 
rotina ZANLYT que pertence ao «package» IMSL (Inter- 
national Mathematical & Statistical Libraries), existente 
no computador DEC1O do LNEC. No caso de (1.7), é 
utilizada a rotina PRBM, já atrás referida. Os polos do 
sistema em anel fechado assim determinados, são repre- 
sentados no diagrama do LGR por um pequeno triângulo, 
caso estejam dentro da área de interesse seleccionada, 


1.4 — Exemplos de aplicação 


1.4.1 — Exemplo | 


Na Fig. 1.4 está representado o lugar geométrico 
das raizes obtido para o sistema 


—s 


e | 
o O(S) = 1: H69) = 


G(s) = 
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Note-se a existência de dois ramos que se iniciam Gts) = s+z 
e terminam na fronteira, cuja representação não seria “ e(s+1) (s+3) 
conseguida se não existisse a pesquisa de fronteira. 


Gls)= 1yHtBj=1 


para Z= 1.5 1.02, 1.01. 


1.4.2 — Exemplo d Na Fig. 1.5 apresenta-se o diálogo realizado entre 


Neste exemplo determinou-se o LGR do sistema o utilizador (caracteres sublinhados) e o programa, no 
caso em que z= 1.5, representando-se na Fig. 1.6 o 
LGR correspondente. 


Nica 
sa 
Eras 
Tindga 
“dna 


CO do E 


14, 


| 
, 
' 
' 
à. 
a 
niva E Eds | 
RR 
Erdra a 


LUG. GEOM. RATZES TIA 


LUG. GEOM.RAIZES 
FIG: ds FIG. 1.5 


NOVO LHCREMENTO od, da 0S* 


SUPERVIGOR UNTY 


VALORES EXTREMOS SELECCIONADOS: 
ndo Us DDOQ0OE HOO 
E ii Liam Sli «EE E deito 
RE nAS= E ÔQONDE-0D 
TAH. 00000E 00 
Epi ro Do Pi DEOMETRIÇO AS RALZES 
UE MONTE TER-LO (& DU NIT = 


PRETENDE RESULTANDO PápioD SISTEMA 
EM AMEI ABERTO? 5 MOU KM o 


NOV SHEN==4 
to METAL pio RaninMal? 
E NOVO YHAX CHMINo=YMAXI=O 


Gogo RALTONAL EM PTEMETERO NO DIS) 
VALORES EXTREMOS SELEBETONADI$ 
LIRC ME PLEHELRONF! INT SMA Da DOE pÓ 
Ensga MINO don GE+HOU 
MM De LDO GE HO 


A DEFTINTU GEO RACLDNAL? ENTHS==2 0 00009E+09 


|= 


PRE TENTE GQ S)c17 6 


Mie PONTOS CilCUlATiDS RE 
do DEFINTO Més RACIQNoOLT HH 
e: PRACA ni DE Li DE BMORT, RELATIVO tão CM HIP ER 


PRETENDE Higor go 


FETA= 5 
das 


REFERENCTA (ATE 5 CARACTERES! T7/SA 

Titan RE LINHA DE AMORT,. KELATIVO Co Qi NI? N 
TRACADO DO LUGAR GEOMETETEO 
DAG RAIZES C'RODT LÓCUS* 3 LOGALTERCAD DE POLOS PARA GANHO 0 €5 OU N9P E 


PUNTOS SINGULARES (REAL Ina.) MH = a 


-0. 1834E +01 O FÓBLE+OO Ro 
O 4 SEDDESDO Ds OUODE OO 


, E PA Pere LOCALTIAÇÃO DOS POLOS PAR mainHio Gs DO000ODE HO! 
DIMENSÕES EXTERIORES DO DESENHO | T [3 Ê POLOS Pai GANHO: S.0000 Õ 


Es e es REL [HA 
E a A AR IUBSAGELOO LS BZOZEAVO 
Hi Tt=14 14 1. 088855409) =. EPTÓTEHOO 
1 BáSTLEADO Us DOQODE TOO 
INCREMENTO CALCULADO! 3, 750E-02 LONA IZACÇÃO DE POLOS FáRA GANHO 0 C& DU N9T MH 


QUER MODIFICAR cs OU MNS* s 
E AGORA ELF! 
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Nas Fig. 1.7 e 1.8 mostra-se o LGR para 2=1.02 e 
2=1.01, respectivamente, 


Note-se as alterações no diagrama quando se modi- 
fica o afastamento entre o zero real e o polo real p= =1. 


LUG. GEOM. RAIZES 


Fig. 1.7 


LUG. GEOM, RAIZES 


FIG. 1.8 


2 — MÉTODO DE NYQUIST 


Este bloco permite ao utilizador obter o diagrama 
de Nyquist [2, 3] de um sistema, conhecida a sua função 
de transferência na forma racional, ou seja, conhecendo 
os valores dos coeficientes dos polinómios do numerador 
e do denominador bem como o valor respeitante ao atraso 
no tempo, caso exista, 

O traçado do diagrama de Nyquist faz-se, fundamen- 
talmente, em duas fases distintas: |) determinação ou 
cálculo dos pontos do diagrama, e ii) traçado propria- 
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mente dito (parte de representação gráfica). Como adiante 
se verá, o algoritmo para a determinação dos pontos a 
representar tem como objectivo evitar zonas de grande 
concentração bem como zonas com diminuto número de 
pontos, o que é conseguido, essencialmente, com base 
no critério de escolha das frequências associadas a cada 
um dos referidos pontos, 


2.1 — Determinação dos pontos a representar 


Quando esta técnica é seleccionada, o utilizador é 
solicitado a fornecer os seguintes dadns: 


— Dimensões exteriores (em cm) do desenho corres- 
pondente ao diagrama; 

— Valores das frequências angulares máxima e mi- 
nima (em rad/8), Wax € Win entre as quais se 
pretende conhecer o diagrama, 


O número máximo de pontos a representar, N, é 
fixado internamente, o mesmo acontecendo com as dis- 
tâncias máxima e minima (em cm), dmax € Amin que se 
admitem entre dois pontos consecutivos. 

Dado que os valores associados a estas distâncias 
limite são estabelecidos em centimetros e que as coorde- 
nadas dos pontos são calculadas em função de unidades 
de G(jw), torna-se necessário conhecer os factores de 
escala (centimetros por unidades de cálculo) relativos aos 
dois eixos para que o critério das distâncias possa ser 
aplicado, Nestas condições, o procedimento preliminar 
adoptado consiste em calcular um número reduzido de 
pontos (cerca de 5% de NJ), distribuidos entre Wi € 
max. é determinar, em seguida, os valores máximo e 
minimo das coordenadas x e y respectivas; com estes 
valores e com o conhecimento, já adquirido, do compri- 
mento total (em cm) dos dois eixos, aqueles factores de 
escala são então estabelecidos, 

Os pontos referentes às frequências limite são sem- 
pre representados. Os valores das frequências de refe- 
rência para o cálculo dos pontos intermédios, isto é, para 
os quais se verifica Win WU <tWa, São tais que 
constituem uma progressão geométrica de N termos em 
que o primeiro e o último são Wi, € WUmaxw respectiva- 
mente, Ter-se-à então que 


(1); > (95.1) (2.1) 


em que a razão, r, é dada por 


In (max Win) 
Ff = exp | MES ÁS rr pe og e a | | 
7 | (2.2) 


O algoritmo respeitante à determinação dos pontos a 
representar, decorre da maneira que a seguir se descreve. 
Para uma determinada frequência tu, são calculadas as 
coordenadas do ponto correspondente; em seguida, atra- 
vês dos factores de escala, é calculada a distância, d(tu;), 
entre este ponto e o último ponto representado. Esta dis- 
tância é comparara com os valores din € dmax podendo 
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êntão verificar-se a ocorrência de um dos seguintes três 
casos: 


1) Amin < d(to,) “ Graz 
2) d(t,) < din 
3) d(ty;) > daax 


No primeiro caso, o ponto calculado é admitido como 
ponto a representar, dado que satisfaz o critério da dis- 
tância, e as suas coordenadas são armazenadas numa 
variável indexada adequada, 


No segundo caso, d(ty) < din, o ponto calculado 
não é representado e passa-se então a outra frequência, 
obtida através de (2.1), repetindo-se o procedimento atrás 
exposto. 

Finalmente, no terceiro caso, d(t);) > das O ponto 
calculado também não é representado e inicia-se um 
processo iterativo de modo a encontrar-se uma frequên- 
cia to; tal que tw; , to; <L to, paraa qual d(w;) satis- 
faça o critério de distância. As frequências ww; são obtidas 
através da expressão 


1/2 
Wj = (14.01) (2.3) 


tendo-se em conta que para j=1 se terá = ty, 
O ponto correspondente à frequência tw, que satisfaça o 
critério de distância é então admitido como ponto a re- 
presentar, O procedimento é então continuado sendo, 
no entanto, de referir que após este processo iterativo 
o primeiro ponto a ser submetido âquele critério é o 
anteriormente calculado e respeitante à frequência tw, 
Por outras palavras, neste terceiro caso, o algoritmo 
permite a representação de tantos pontos, corresponden- 
tes a frequências no intervalo (to, ,, tu;)), quantos os ne- 
cessários para que o critério de distância seja satisfeito. 


2.2 — Representação gráfica 


Os valores das coordenadas dos pontos a represen- 
tar são transmitidos a esta fase na qual o diagrama de 
Nyquist é traçado, utilizando-se para tal um algoritmo 
que não apresenta quaisquer caracteristicas especiais. 

Nesta fase, o utilizador pode ou não optar pela gra- 
duação do diagrama com valores da frequência angular. 
Esta graduação é feita pela colocação de marcas sobre 
o diagrama, com indicação do valor da frequência corres- 
pondente. Para se efectuar esta graduação utilizam-se os 
valores armazenados em «arrays» apropriados que con- 
têm os valores da frequência e das coordenadas do 
ponto correspondente, valores estes obtidos durante a 
determinação dos pontos a representar, Os valores cor- 
respondentes a Win & Wax São automaticamente consi- 
derados para este efeito; para além destes, os pontos 
utilizados para a graduação são aqueles cuja ordem de 
obtenção é um múltiplo de um número, internamente 
fixado, que é cerca de 109% do número máximo de pon- 
tos a representar, 
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2.3 — Exemplos de aplicação 


2.3.1 — Exemplo | 


Neste exemplo determinou-se o diagrama de Nyquist 


do sistema 


G(s) = 


EX 


é CO. Got; HO) 1 


para a =0O(Fig 21) ea =1 (Fig. 2.2). 


DIAGRAMA NYQUIST 


= 
— 


FIG. 2.1 


Note-se a modificação do diagrama resultante da 
introdução de atraso no sistema, 


DIAGRAMA NTQUIST 


FIG. 2.2 


2.3.2 — Exemplo Il 


Nas Fig. 2.3 e 2.4 mostra-se o diagrama de Nyquist 
obtido para o sistema 


G(s) = 


80 


s3 + s? + 100s 


;C(s)=1:H(s8)=1 
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e parte do diálogo realizado entre o utilizador (caracteres 
sublinhados) e o programa, 


DIAGRAMA NYQUIST 


FIG. 2.3 
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FIG. 2.4 


3 — METODO DE BODE 


Este bloco permite ao utilizador obter o diagrama 
de Bode [2, 3] de um sistema, conhecida a sua função 
de transferência, No nosso caso, a implementação deste 
método implica o conhecimento não só do valor do atraso 
no tempo mas também dos coeficientes dos polinômios 
do numerador e do denominador daquela função na sua 
forma racional, 


Tal como para as técnicas já atrás referidas, o algo- 
ritmo adoptado para a determinação dos pontos a repre- 
sentar no diagama de Bode tem em especial consideração 
o objectivo que consiste em evitar a desigual distribuição 
daqueles pontos. Ver-se-á adiante que este objectivo é 
conseguido através de uma escolha adequada das fre- 
quências correspondentes a esses pontos, 


3.1 — Determinação dos pontos a representar 


O diagrama de Bode é constitudo por duas curvas 
que estabelecem o andamento do módulo e do ângulo 
de fase de G(jw)* em função da frequência angular 
tu(rad/s), Para o módulo e para a frequência usam-se 
escalas logaritmicas (na base decimal). Relativamente ao 
módulo, este expressa-se em unidades designadas por 
decibeis (dB), de tal modo que 


| G(juo) | [dB] = 20 log | G(jto) | (3.1) 


A estrutura do programa obriga a que o dominio da 
frequência seja constituido por um número inteiro de 
décadas de frequência, Assim, quando esta técnica é 
seleccionada, o utilizador é solicitado a fornecer o número 
de décadas, N,, que pretende; para além deste dado, 
deve ainda indicar o número de pontos a representar, N, 
e o valor da frequência angular inicial, tw». O valor da 
frequência inicial deve ser uma potência positiva ou ne- 
gativa de 10; caso o não seja, o algoritmo escolhe, auto- 
maticamente, o valor imediatamente inferior que satisfaça 
aquela condição, dando deste facto indicação ao utili- 
zador. 

Tendo em atenção o objectivo atrás mencionado, as 
frequências de cálculo ty, são determinadas pela expres- 


são seguinte 
Na 0—1) — 
| (N—1) | 


td, = 405.10 (= 12, 7 N) (3.2) 

Para cada uma desta N frequências são calculados 
os valores do módulo de G(jtw), em dB, e da fase, em 
graus, valores estes armazenados em variáveis indexadas 
adequadas. 


3.2 — Representação gráfica 


O algoritmo utilizado para a representação gráfica 
dos valores obtidos como atrás se indicou não apresenta 
quaisquer caracteristicas especiais. Há, no entanto, que 
referir que nesta fase o utilizador deve indicar as dimen- 
sões exteriores (em cm) desejadas para o desenho do 
diagrama de Bode. São também indicados ao utilizador 
os valores máximo e minimo seleccionados para o módulo 
e para a fase, os quais podem ser alterados, caso tal 
seja pretendido. 


3.3 — Exemplos de aplicação 
3.3.1 — Exemplo | 


Neste exemplo determinou-se o diagrama de Bode 
do sistema 


100 —8 


G(g)-=-———————— e - 
(8) sº + 1.5s + 100 


cs) =1; H(g=1 


(*) Para simplicidade da exposição, designamos por Gí(jw) quer a função de transferência em anel aberto (produto CGH) quer a 


função de transferência em anel fechado. 
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para 4x =1 (Fig. 31) e «=0 (Fig. 3.2). 


Note-se a alteração na fase do sistema resultante da 


introdução de atraso. 


Bt 


DIAGRAMA BODE 


FIG, 3.1 


á! 


DIAGRAMA ECDE 


FIG. 3.2 


3.3.2 — Exemplo dl 


Nas figuras 3.3 e 3.4 mostra-se o diagrama de Bode 


obtido para o sistema 


500 s“ + 500 s +40 000 


( G|Dl —=s == — = = 
(s) s!+11.5 s!4+115 52410005 


e parte do diálogo realizado entre o utilizador (caracteres 


sublinhados) e o programa, 
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: C(s)=1; H(g)=1 


FIG. 3,3 


REFERENCIA (ATE 5 CARACTERES): BS/5A 
NR. PONTOS DE CALCULOCS=1000)*? NPT=500 


NkK. DE DELADAS DE FRED. (>070=5)7 NDEC=-? 


FREQ, ENECIALCMULT. OU SUBMULTO ME 157 MICRADAS)-Laa 
TOMOU -SE PARA FREG. INTCIAL WILRAD/G)= 1, 0E+400 
TRACADO DO DIAGRAMA DE BODE 


DIMENSDES EXTERIORES Dil DESENHOS 
EA=< X io à B E é=I9 4EMII 
Mof=1i72 19 


JALURES EXTREMOS CALCULADOS 
REL «AMPI.I TUDNES (hECTIREISO: 
Hink S 1FREOI 
AIN=-2L0SEHOL 

Foo CORALGOs 

(A UI GE FO 

Hind= 1].7"40E+OZ 

QUER MODIFICAR (5 OU MT E 


REL. AMPLITUDES (DECTRETSS! 
NOVIS MAX=40 
NOVO MINc=2a 


MEO MIM DEVE GER So MINI 
MO MT A) 


FASES CORBIS 


HNUVO MAX=Q 
IVO MEN=-=1H0 


ADORA? SUP 


FIG. 3.4 


4 — RESPOSTA NO TEMPO DO SISTEMA 


4.1 — Introdução 


Este bloco permite determinar a resposta no tempo 
em anel aberto ou fechado de qualquer sistema, conhe- 
cidas as funções de transferência dos subsistemas cons- 
tituintes, nomeadamente, do sistema controlado, do con- 
trolador e do sistema de retroacção, 
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A sequência de cálculo adoptada está indicada, de 


uma forma muito simplificada, no fluxograma da Fig. 4.1. 


RT EMP O 


(RESPOSTA 
NO TEMPO! 


“Entrada 
dados 


FT — " MBECÇO 


Resposta no tempo] 


do sistema 


; = = 
Representação 
gratica 


MNEMONIC A 
“ERRADA 


Mensagem | | 


Mensagem 


explicativa crio 


FIG. 4.1 


A simulação do sistema é realizada com base na sua 
representação em Espaço de Estados, sendo esta deter- 
minada pela subrotina RMIN, já definida no Bloco de 
Manipulação de Dados ([1]. 

A integração das equações diferenciais da represen- 
tação na forma de Espaço de Estados do sistema, é 
realizada pela rotina RKM, que utiliza o método de Runge- 
-Kutta-Merson [6]. 


4.2 — Resposta no tempo do sistema 


Para determinar a resposta no tempo do sistema, é 
necessário definir um conjunto de parâmetros adicionais, 
indispensáveis à integração da sua representação em 
Espaço de Estados. 
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Estes parâmetros, introduzidos pela sub-rotina 
DSIMUL, chamada na sub-rotina RKM, são os seguintes: 


q) — tempo inicial 

“TE — tempo final 

- DT —- incremento no tempo 
*H — passo de integração 
* FORM — forma da entrada 


São ainda estabelecidos pelo programa o erro abso- 
luto máximo por segundo (Emax = 0.00001) e o valor ini- 
cial para as variáveis de estado (x, = 0, VY;) dando-se 
a possibilidade ao utilizador de alterar estes valores, 
caso o pretenda. 

A representação gráfica da resposta do sistema no 
tempo — Y(t) — será construida com base nos valores 
calculados para a saida nos instantes 


ti; = t..; + DT 
impondo-se como restrições de programação 


105 2H<£DT 
NPT < 1000 
em que 


NPT —n.º de pontos de cálculo (NPT=1+(TF—TI)/DT) 


4.2.1 — Tipo de entrada 


Admitem-se as seguintes possibilidades para a forma 
de variação no tempo da entrada; 


* DEG — entrada em degrau 

* RAM — entrada em rampa 

* PAR — entrada parabólica 

* SIN — entrada sinusoidal 

* VAR — entrada variável, definida numa tabela. 


4.2.1.1 — Entrada em degrau 


U; para t 2 TI 
T) = 
u(t) o para t < TI 
em que 
u;— constante real 
4.2.1.2 — Entrada em rampa 
Ro, (t-TI) para t > Tl 
u(t) = | 
O para t< TI 


em que 
R,— declive da função rampa (constante real positiva) 
4.2.1.3 — Entrada parabólica 


Ro (t-TI)? para t > Tl 
ult) = 
O para t < TI 
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em que 


AR, — constante real positiva 


4.2.1.4 — Entrada sinusoidal 


u; + À sin[w(t-TI)] para t > TI 
t) = 
u(t) O para t < TI 


em que 
u; — constante real 


A — amplitude 
tó — frequência angular 


tigd4.? Entrada em degrau 


Higd4A4 Entrada parabolica 


Os valores de u(t) correspondentes a instantes não 
definidos na tabela são obtidos por interpolação linear, 
tendo-se: 


(U, — Up). (t— Tua) 


u(t) = 


+ U,., para Ty <t<T, 


em que 
U, — valor tomado pela entrada no instante T, 


4.2.2 — Integração de equações diferenciais pelo mé- 
todo de Runge-Kutta-Merson [6] 


O método utilizado na simulação de sistemas lineares 
univariáveis é o de Runge-Kutta-Merson de 5.º ordem, 


“ig 1.3] Entrada em rampa 


fig 4) Entrada sinusoidal 


FIGS. 4.2, 4.3, 44 0 4.5 


4,2.1.5 —- Entrada variável 


Os valores tomados pela entrada e correspondentes 
instantes no tempo são definidos por uma tabela, sendo 
os pontos dados por ordem crescente do tempo, e obe- 
decendo às seguintes condições: 


 TNOCÇT 
* T(NP) > TF 
* 24 NP<Z100 


em que 


NP — n.º de pontos da tabela, 


que apresenta como vantagem fundamental a possibili- 
dade de se controlar o passo de integração com base 
numa estimativa do erro cometido, 

Consideremos um sistema de equações diferenciais 
ordinárias 


dx 
dt 


=g (x, u) 


A solução deste sistema de equações tem a forma 
geral 


5 
x (t+Hh)=x(0+H wk; (4.1) 
j=1 
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ém qué 
h — passo de Integração 


k; — incrementos calculados em diferentes pontos 
w; — pesos dos incrementos 


Admitindo para os pesos dos incrementos, os valores 


Wj = 1/6 
io = W3 = O 
wa = 4/6 


a expressão (4.1) assume a forma 


x (t+h) = X (O) + (k; + dks + ks) /6 


sendo 
k; =h. g (x + U), 
ks ae UP q (x + k; [3 U)ieh/3 
k=h.g (x+ (k+k)/6., Una 
k=h.g (x+ (kj+3k3)/8., U)un/> 


ks=h.g (x+ (kj—3k+4k4)/2., U)n 


Designando por E = [£;] a estimativa do erro come- 
tido durante um passo de integração, tem-se 


e=(2k—9 k+8 k—ks)/30 (4.2) 


Designando por E, = [E] a normalização do erro 
relativamente ao passo de integração, tem-se 


E; 
; 


A variação do passo de integração resulta da compa- 
ração de £,; e Enax de acordo com o seguinte quadro: 


Mantêm-se o passo de in- 
tegração constante 


N= 28 
Considera-se válido o pon- 
to calculado 


H' = Hh2 
Não se considera válido o 
ponto calculado 


3, - Eni > Eni 


4.3 — Exemplos de aplicação 


Nas figuras 4.6, 4.7, 4.8, 4.9, 4.10 mostra-se a res- 
posta no tempo obtida, em anel aberto, para o sistema 


156-25 


(s) sº + 25 s + 156.25 
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C(s) = 1 
H(s) = 1 


Nota: A tracejado foi representado, pelo programa, 
o sinal de entrada do sistema, 


FIG. 4.6 


Resposta do sistema so degrau unitário - 


RESPOSTA NO TEMPO 


à aa —— — Te—e— 


FIG. 4.7 


Resposta do sistema ao sinal rampa unidade 


RESPOSTA NO TENPÓ 


FIG. 4.8 


Resposta do sistema ao sinal parábola unidade 
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RESPOSTA NO TEMPO 


FIG. 4.5 


Resposta do sistema ao sinal sinusoidal de ampii- 
tude unitária e frequência igual à frequência própria 
do sistema 


RESPOSTA NO TEMPO 


FIG. 4,10 


Resposta do sistema à onda quadrada de amplitude 
unitária e periodo igual a i0s 


Na figura 410 o sinal de entrada (onda quadrada) 
foi gerado por pontos utilizando a facilidade criada pela 
entrada variável (VAR), 

Na figura 4.11 apresenta-se uma parte do diálogo 
entre o utilizador (caracteres sublinhados) e o programa, 


realizado na obtenção da resposta do sistema à entrada 
sinusoidal (Fig. 4.9) 


GEFERENCIA CATES 5 CARACTERES)! KTL 


à RESPOSTA NO TEMPO SERA OBTIDA PARA OS INSTANTES T(N)+ TAIS QUE! 
TLHI=TOM=LI+0DT 
SAO DADOE é 


Li INCREMENTO NÚ TEMPO 
TI TEMPO INICIAL 
TF TEMPO FINAL 


TEMPO INICIAL =Q, 
TEMPO FINAL =g 

INCREMENTO KO TEMPO =.01 
PASSO DE INTEGRACAO =,005 


[ERRO ABSOLUTO MAXIMO =1,E-1) 
FRETENDE MODIFICAR" A 


ERRO ABSOLUTO HAXIMO= +01 


CUARIAVEIS DE ESTADO INICIAIS NULASI 
PRETENDE MODIFICAR?H 


TIPO DE ENTRADA: SIM 


ENTRADA SINLUSOIDAL 


U=UL+HAMP. SINCWLT=TI AI 


DADOS: 

IT ENTRADA INICIAL 

u FREQUENCIA ANGULAR 
AMP AMPLITUDE 
UINICIAL= 

AMPLITUDE =] 


FREQUENCIA ANGULAR=12,5 


INICIADA INTEGRAÇÃO (R=k=M) 


FIM DA INTEGRACÃO 


TENCADO DA RESPOSTA NO TEMPO 


PIMENSCES EXTERIDRES DO DESENHO 
CB=<. X <=30 à B=€< T <=30 (CMI 


XeY=12 2 


E AGORA? SUP 
FIG. 4.11 
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Tecnologia da metanogénese 


(Revisão descritiva) 


RESUMO 


Neste artigo caracterizam-se as variáveis actual- 
mente consideradas como fundamentais no processo 
metanogenético, descrevem-se as diferentes arquitectu- 
ras mais correntes nos digestores rurais, industriais e 
laboratoriais, sem que se faça qualquer análise quanti- 
tativa das relações matéria-prima — produtos. No en- 
tanto procura-se tanto quanto possível apresentar uma 
bibliografia exaustiva do estado actual da tecnologia 
de metanogénese. 


1 — INTRODUÇÃO 


Nos últimos anos o «aparecimento da crise ener- 
gética» levou à necessidade do estudo das fontes de 
energia, e assim ressurgiu a METANOGENESE, isto é, 
o estudo da decomposição por via anaeróbia dos com- 
postos orgânicos [1, 2, 3, 4, 5], com produção simul- 
tânea de um gás combustível: o BIOGAS (constituído 
por aproximadamente 60 % de metano). 


O esquema da figura 1 sintetiza o processo da 
metanogénese. As matérias-primas (compostos orgâni- 
cos animais ou vegetais) fornecidas ao digestor (reac- 
tor biológico) são aí retidas durante algum tempo 
(tempo de residência), decompostas e transformadas 
nos produtos finais pela biomassa presente (população 
bacteriana). Esses produtos finais são: 


— um gás combustivel (o biogás); e, 
— um resíduo liquido-sólido, que se utiliza como 
fertilizante. 


MARIA CELESTE ANTUNES () 
SILVESTRE DIAS ANTUNES (")(*) 


ABSTRACT 


In this article, variables nowadays considered as 
fundamental in the anaerobic digestion process are cha- 
racterized, and the different architectures of digesters 
used in rural, industrial and laboratory applications are 
described. Although, if quantitative analysis of the 
transformation: prime matter — products is not consi- 
dered, an exhaustive bibliography on the present state 
of the anaerobic digestion technology has been tried. 


MATERIAS 


PRODUTOS 
PRIMAS INAIS 


| METANOGENESE 


[Compostos Ê ertihzante em 
orgânicos em | suspensão aquosa, 
suspensão aquosa | e biogás ) 


DISESTOR 
(reactor biológico, Contendo 
o bromassa ) 


FIG. 1 


O processo da metanogénese 


A metanogénese com o objectivo exclusivo da 
produção energética de metano poderá não ser um 
processo economicamente viável. No entanto, a utiliza- 
ção da metanogénese é essencial para a modificação 
do bem-estar das populações que vivem em zonas não 
servidas pelas redes de distribuição de energia, permi- 
tindo assim: 


() — Licenciada em Farmácia (U. L.), estudante de pós-graduação na U. C. L. (Bélgica), 
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(*) — Morada actual: 


Université Catholique de Louvain 
Unitá Auto 

Bátiment Maxwell, Place du Levant 
1348 — LOUVAIN-LA-NEUVE 
Bélgica 


Manuscrito recebido para publicação em 4/1/B0. 
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-=— a produção de um gás combustível utilizável 
no aquecimento doméstico, ou na produção de energia 
eléctrica (recorrendo à utilização de motores de com- 
bustão interna em instalação fixa, do tipo do grupo 
TOTEM da fábrica FIAT [12]); 


-—— a produção de adubos não-sintéticos; 


— o tratamento de resíduos orgânicos (animais 
bu vegetais) focos de poluição; e, 


-— | fomento e respectivo desenvolvimento de 
uma produção animal e agricola (fontes da matéria- 
-prima para a produção do biogás e dos adubos não- 
-sintéticos). 

Aproveitamentos quer colectivos [7], quer indi- 
viduais [6] têm sido efectuados e desenvolvidos com 
bons resultados em países onde ainda não existem 
redes de distribuição de energia, tais como a China 
e a India. Nesses casos o digestor ideal do tipo faça-o 
você mesmo, é construído com materiais (localmente) 
acessíveis, segundo técnicas que os governos põem à 
disposição das populações. 

Em locais melhor servidos energeticamente (aglo- 
merados populacionais, industriais,...) é imperioso o 
tratamento dos focos de poluição [8, 9] (tratamento 
da lixos, esgotos, resíduos industriais,...) e a metano- 
génese permite então aproveitar o biogás produzido 
como factor de recuperação dos investimentos efec- 
tuados nas instalações de tratamento. 


Actualmente começam a surgir estudos [10, 11] 
de utilização da biomassa na produção de combusti- 
veis líquidos, obtidos por diálise da fase líquida con- 
tida nos digestores. 


Na figura 2, sintetiza-se esquematicamente o pro- 
cesso biológico da metanogénese (um dos esquemas 
actualmente proposto e em discussão pelos biólogos, 
ver referências [1] a [5]). 


COMPOSTOS 


ACIDOS 


MATERIA BRiMa = ORGÂNICOS ORGÃ did 
Compostos orgânicos) |” | (acel , COMPLEXOS 
[Compostos orgânicos) | (aeetico, | fceluloses 


SPOPI ÓCIO, 
F hemuce lubires 


q - Ea 
ips LF 


produtos qua dos 


HO 


FIG. 2 


Transformação biológica da matéria-prima 


Nos parágrafos seguintes vamos ver como se ca- 
racterizam actualmente as matérias-primas existentes, 
estudar as características dos produtos da metanogé- 
nese e comparar a arquitectura dos diferentes tipos de 
digestores existentes. 
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2 — AS MATÉRIAS-PRIMAS 


As matérias-primas de origem animal ou de origem 
vegetal após trituração e mistura adequada, ou pré-tra- 
tamento térmico, aeróbio, químico ou sónico, são intro- 
duzidas no digestor em suspensão aquosa (a percenta- 
gem de água depende essencialmente do tipo de 
digestor). 

Inicialmente utilizavam-se apenas como matérias- 
“primas os estrumes resultantes das explorações agri- 
colas-animais. Actualmente faz-se também o aproveita- 
mento dos resíduos industriais não só com o objectivo 
da produção do biogás, mas essencialmente com o fim 
de tratar os efluentes industriais. 

A caracterização das matérias-primas é complexa, 
e faz-se quantitativa e qualitativamente. 


214 — Caracterização das 
quantitativo 


matérias-primas: aspecto 


Ao alimentar um digestor define-se um caudal F 
de entrada da matéria-prima (volume por unidade de 
tempo) e um tempo Q de residência dessa matéria no 
digestor (que deve ser superior ao tempo necessário 
para a transformação da matéria). Se for V, o volume 
da fase líquida do digestor, esse tempo será: 


0= (2.1.1) 


Este tempo de residência da matéria no digestor, 
é condicionado pela velocidade de transformação da 
matéria, e esta depende da quantidade e tipo da bio- 
massa presente, do tipo de digestor e da qualidade da 
matéria-prima. 

No aspecto quantitativo podemos considerar três 
tipos de alimentação das matérias-primas: descontínua, 
intermitente e contínua, 

Em digestores descontínuos (tipo batch), as ma- 
térias-primas são carregadas em simultâneo com a bio- 
massa (resíduo parcial de uma digestão anterior) e a 
transformação processa-se sem alimentação posterior, 
até que deixe de haver produção do biogás, ou até que 
esta seja minima, após O que se procede a um novo 
carregamento do digestor. 

A alimentação intermitente (digestores contínuos) 
é vantajosa nos casos em que a obtenção da matéria- 
-prima também é intermitente, evitando-se nesse caso 
a necessidade de uma armazenagem inicial, É o que 
se verifica nas instalações agrícolas-animais onde os 
resíduos da exploração são recolhidos diariamente, e 
diariamente introduzidos no digestor. 

A alimentação contínua (digestores contínuos) é 
importante no caso do tratamento de efluentes prove- 
nientes de indústrias de produção contínua. 


2.2 — Caracterização das matérias-primas: aspecto qua- 
litativo 


A determinação da qualidade de uma matéria- 
“prima é difícil. Vários métodos têm sido utilizados, 
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dependendo essencialmente da natureza dos substractos 
consumidos pela biomassa presente no digestor. Veja- 
mos cada um dos métodos: 


2.2.1 — Análise química ponderal dos elementos 
químicos mais importantes: carbono, hidrogénio, oxi- 
génio, azoto, fósforo e enxofre presentes na matéria- 
-prima. Determina-se deste modo uma fórmula química 
aproximada C.H,O.N,;P.S, com base na qual poderão 
determinar-se [13, 14] os produtos químicos das reac- 
ções hipoteticamente realizadas no digestor. 

Atendendo a que o carbono e o azoto são elemen- 
tos fundamentais, utiliza-se correntemente como coefi- 
ciente de qualidade de uma matéria-prima a proporção; 
carbono/azoto (C/N) indicando-se como valor ideal 
para a digestão o valor 30/1. Ainda que o valor 30/1 
seja o ideal, a matéria-prima poderá apresentar uma 
digestão mais ou menos difícil (demorada) consoante 
a forma sob a qual se apresentam os diferentes ele- 
mentos. 


2.2.2 — Análise física ponderal por aquecimento 
prolongado a temperaturas constantes: a análise física 
ponderal após aquecimento prolongado a temperatura 
constante permite a determinação das matérias que se 
volatilizam dentro de determinadas gamas de tempe- 
ratura. Assim utilizam-se correntemente [6] as seguin- 
tes definições: 


— SÓLIDOS TOTAIS (ST) — composição ponde- 
ral da matéria-prima em sólidos residuais após 
aquecimento prolongado (aproximadamente 4 
horas) a 105ºC (secagem total); 

— SÓLIDOS VOLATEIS (SV)] — composição pon- 
deral da matéria-prima em sólidos que se vola- 
tilizam após secagem total (isto é, após deter- 
minação dos sólidos totais) por aquecimento 
prolongado (aproximadamente 8 horas) a 
450 “C; 

— CINZAS (C,.) —- composição ponderal da ma- 
tária-prima em sólidos residuais após aqueci- 
mento prolongado (aproximadamente 3 horas) 
a 1000º C; e, 

— CARBONATO DE CÁLCIO (C,CO,) — compo- 
sição ponderal da matéria-prima em sólidos 
volatilizados entre os 450º C e os 1000ºC, 
expressa em carbonato de cálcio: 


100 


C,CO, = (ST-SV-C,) . “RE [Asas 


A escolha da técnica de determinação da quali- 
dade é discutível, na medida em que, a análise química 
pondera! é extremamente complexa e conduz a resul- 
tados pouco reais. A análise física ponderal parece ser 
mais real dado considerarem-se actualmente os sólidos 
voláteis [15, 16] como fontes mais prováveis do me- 
tano produzido pela biomassa. 

No caso dos efluentes industriais utiliza-se corren- 
temente como característica de qualidade da matéria- 
“prima um dos parâmetros: demanda quimica (COD) 
ou bioquimica (BOD) do oxigénio necessário para a 
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oxidação total das matérias transportadas pelos referi- 
dos efluentes. 

No aspecto de qualidade determina-se ainda o pH 
da matéria-prima, factor considerado como um dos 
mais importantes para a sobrevivência da biomassa 
(meio ambiente da biomassa), e/ou para que se não 
produza inibição da sua actividade. 


Como qualidade da matéria-prima deve ainda con- 
siderar-se a presença de biomassa activa. Em parti- 
cular, nas matérias de origem animal (como no caso 
dos excrementos dos ruminantes). Actualmente uma 
só técnica [17] permite caracterizar a fracção meta- 
nogenética dessa biomassa: a análise do coenzima Fog: 
que dá uma indicação aproximada da biomassa meta- 
nogénica activa presente na matéria analisada. 


3 — BIOMASSA 


Designa-se por biomassa o conjunto das popula- 
ções bacterianas activas que permitem a decomposição 
da matéria-prima com produção simultânea de biogás. 
No esquema da figura 2 apresentou-se o conjunto das 
diferentes fases consideradas no tratamento anaeróbio 
dos compostos orgânicos. 

A classificação, caracterização e funções da bio- 
massa, tal como a sequência efectiva de decomposição 
da matéria, não é ainda bem conhecida. Várias classi- 
ficações têm sido apresentadas pelos biologistas (ver 
referências [1] a [5]). E difícil fazer uma análise 
termodinâmica das reacções complexas desenvolvidas 
durante a metanogénese, pois que a cinética de for- 
mação de certos produtos intermédios (tais como áci- 
dos orgânicos em C,, C;,...) inibidores da acção me- 
tanogenética é difícil de prever [18]. 

A actividade da biomassa presente (entenda-se por 
actividade biológica a formação de nova biomassa com 
produção de produtos secundários no digestor) depen- 
de, das condições do meio ambiente da biomassa: 


— do pH, e este do poder tampão do meio; 

—— da temperatura ambiente; e, 

— da presença de produtos químicos (tóxicos, ini- 
bidores e/ou estimulantes). 


A temperatura é talvez o factor ambiente funda- 
mental dado que o tipo da biomassa metanogénica 
depende desta. Consideram-se [19] normalmente três 
tipos de biomassa metanogénica de características di- 
ferentes consoante a temperatura a que desenvolvem 
a sua actividade máxima. Assim existem: 


— uma biomassa criófila — cuja actividade é má- 
xima a temperaturas de 12 a 18º €; 

—— uma biomassa mesófila — cuja actividade é 
máxima a temperaturas de 25 a 40º C; e, 

—— uma biomassa termófila — cuja actividade é 
máxima a temperaturas de 55 a 65ºC. 


As biomassas mais utilizadas são de natureza me- 
sófila (por razões económicas — aquecimento dos di- 
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gestores), ainda que a utilização de biomassas termófi- 
las pudesse ter vantagens, nomeadamente na eliminação 
de organismos patogénicos e na degradação de com- 
postos orgânicos complexos. 


4 — OS DIGESTORES (REACTORES BIOLÓGICOS) 


O digestor ou reactor biológico ideal será aquele 
em que toda a matéria-prima introduzida é transforma- 
da em biogás, e cujos produtos finais sejam um líquido 
isento de matéria-prima e de biomassa. 

A impossibilidade da transformação total da ma- 
téria-prima por via anaeróbia, permite a utilização do 
resíduo sólido-líquido como fertilizante na agricultura, 
quer directamente, quer após tratamento adequado, 
O fertilizante obtido por via anaeróbia tem sobre o 
obtido por via aeróbia, a vantagem de conter uma 
menor percentagem de micro-organismos vivos e de 
conter uma maior percentagem de azoto. 

Existem diversos tipos de digestores, cuja forma 
e tecnologia dependem de diversos factores: 


— a localização e temperatura ambiente, que le- 
vam à necessidade de aquecer ou não os diges- 
tores consoante a natureza da biomassa utili- 
zada; 

— a natureza da matéria-prima utilizada (efluente 
industrial, resíduo agricola-animal,...); 


— modo de alimentação do digestor (continua, 
descontínua ou intermitente); 
— estrutura da matéria no digestor: homogénea 


ou heterogénea ([digestor com ou sem agita- 
ção); e 

—— tipo do digestor: rural ou industrial (eventual- 
mente laboratorial). 


Vamos apresentar a arquitectura de alguns tipos 
de digestores. Assim consideraremos: 


— digestores convencionais: 


e do tipo descontínuo (batch); e, 


do tipo contínuo de alimentação inter- 
mitente, 


— digestores com realimentação de sólidos de- 
cantados; 

— digestores com fases cineticamente separadas 
(acidificação e metanogénese separadas); 

— filtros anaeróbios; 

— digestores do tipo de escoamento vertical 
fup-flow]; e, 

— digestores heterogéneos do tipo de escoamento 
sólido fplug-fowa. 


No quadro | indicam-se as matérias-primas mais 
utilizadas para cada tipo de digestor, e as figuras 3, 4 
e 5 representam convencionalmente alguns dos tipos 
de digestores. 


QUADRO | 


Tipos de digestores e matérias-primas mais utilizadas 


TIPO DE DIGESTOR: 


Digestor descontinuo 


convencional | contínuo 


(= [DDD [DD == os mess... = ss; a 


Digestor com realimentação de sólidos decantados 


Digestores com fases cineticamente separadas . 


Filtro anaeróbio 


Digestores do tipo up-flow 


a 


Digestores do tipo plug-flow... 


TIPO DE MATÉRIA-PRIMA: 


Efluentes 
industriais 


Resíduo de exploração 
agricola-animal 


cH, + CO 


Matérias 
primas 


E 
—= |———— efluente 


É 


DIGESTOR CONVENCIONAL 
FIG. 3 


Digestor convencional homogêneo 


S12 


Naférias 


pri mas 


E Recirculação de solidos | 


FIG. 4 


Digestor convencional com realimentação de sólidos decantados 
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Matérias primas 


FILTRO ANAERÓBIO 
FIG. 5 


Filtro anaeróbio 


4.1 — Digestores convencionais 


Os digestores convencionais são de dois tipos. 

O mais antigo, e mais conhecido é o digestor 
descontinuo (ou de tipo batch) em que a matéria-prima 
é introduzida em conjunto com a biomassa num reser- 
vatório, desenvolvendo-se a produção do biogás durante 
um período de tempo mais ou menos extenso, função 
da dimensão do reservatório (isto é, da quantidade de 
matéria-prima) e da quantidade de biomassa inicial. 

Segundo Lagrange, B. [20], existem em Portugal 
três digestores deste tipo, que funcionam sem proble- 
mas desde há vinte anos. 

O manual da W. P. C. F. [29] apresenta alguns 
desenhos de digestores do tipo descontínuo, em betão 
armado; 


— digestores com cobertura em aço, fixa, porém 
amovível (figura 6); e, 

— digestores com cobertura em aço, móvel (fi- 
gura 7). 


Mas poderão ainda considerar-se o caso de diges- 
tores com cobertura deformável, funcionando esta como 
reservatório para gás (figura 8) [26]. 


FIG. 6 


Digestor em betão armado com cobertura fixa, 
amovivel [29] 


E COMATURA novEi , 
À EMA  — 


FIG. 7 


Digestor em betão armado com cobertura 
móvel [29] 
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fame ch 
| / did RCA. EFLUENTE 
PA E vm» a a É 
DE AunentAÇÃO A as o Ge fá 
DISPPEPIES SEER EANES TS. 
FIG. 8 


Digestor com cobertura deformável, funcionando como reservatório 
para o biogás [25] 


Os mais utilizados, dada a sua grande expansão 
nas pequenas explorações agrícolas-animais são os di- 
gestores convencionais de alimentação intermitente. 

Nas figuras 9 a 23 apresentam-se várias arquitec- 
turas de digestores do tipo convencional, indicando-se 
sempre que possível a sua origem e referências biblio- 
gráficas. 


FIG. 9 


Digestor convencional descontinuo, tipo indiano [22] 


quto PE AtucuiaçÃo 
pa FASE LIQUIDA 


CtTo En MPE E Tivbo 
[lº TEU ME ( Fase sElida) 


paro: Ay] 
AI sonma DE 
ANO |  ErEmEusAÇÃO. 


FIG. 10 


Digestor convencional descontinuo, tipo |. N.R.A. (França) [21] 
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Projecto Industrial de um digestor convencional descontinuo com 


E FIG. 14 
agitação e bombagem em circuito fechado da matéria liquida de 
modo a evitar a formação de uma camada de espuma que possa Esquema da instalação de tratamento de estrumes em digestor 
impedir a saida do biogês, segundo Horton, R. [33] convencional com carga-descarga motorizada [35] 
É À peouás 
E: 
RA | 
E GASONETRO à 
L 
E f aa 
1 
+ 
sala cs =D === = == 
Ai) 


FIG. 12 


Digestor convencional com alimentação Intermitente e armazenamento do gás 
em gasômetro (India) 
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PARA BIDGAS 
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FIG 13 FIG. 15 
Projecto industrial de um digestor do tipo convencional para Forma típica do digestor convencional de alimentação intermi- 
alimentação intermitente ou continua tente, tipo indiano [30] 
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tipo chinês [32] 


FIG. 17 


Forma típica de uma instalação 

de produção de biogás do tipo 

chinês, segundo a FAO Soils 
Bulletin [32] 


FIG. 18 


Exemplo de uma instalação in- 

tegrada tipo chinês [32], para 

tratamento de residuos agrico- 
las-animais e domésticos 
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o 


GASÓMETEO 


à Ta 
ENTRADA prog Shi JD 


SAÍDA HG 
ANA 
| E E E T Tratamento e medida do 
cÂnAMA de FERMENTAÇÃO 


Broa j 
E usas La ana 
FIG. 19 


Digestor convencional com agitação continua ou intermitente, de tipo semelhante ao chinês, segundo 
Horton, R. [33] 


Fi 


A Pd Fl r, 
é EE / 
FIG. 20 


FIG. 21 
Digestor chinês tipo esférico, de pequenas dimensões [31] 


Técnica de construção do digestor chinês, 
tipo esférico [31] 


q—- 
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FIG. 22 
Digestor tipo saco (saco em neoprene) utilizado em Taiwas, segundo Pyle, D.L. [34] 
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FIG. 23 


Digestor convencional de alimentação in- 
termitente, tipo laboratorial (601 de volu- 
me util), segundo Nyns, E, J. el al, [30] 


4.2 — Digestores com fases cineticamente separadas 


Segundo Gosh, S. et al. [36] a digestão em duas 
fases cineticamente separadas (acidogénese e metano- 
génese) apresenta vantagens e desvantagens. Assim: 


—— melhora-se a possibilidade de manter para cada 
grupo de bactérias (acidificantes e metanogé- 
nes) as condições ideais de actividade, ainda 
que actualmente tal separação seja considerada 
prejudicial pois considera-se [1, 2, 3, 4, 5] 
que as referidas bactérias actuam em sim- 
biose; 


— reduz-se o custo do digestor dado que o volume 
total é menor; 


—— aumenta-se a produção de metano: e, 


— reduz-se a quantidade de energia utilizada no 
aquecimento do digestor, pois só a fase meta- 
nogénica deve ser aquecida. 


= 


Apesar de tudo, este processo é mais caro pois 
carece de mais instrumentação e além disso a instala- 
ção é mais difícil, ainda que a separação seja ciné- 
tica. 


Na figura 24 representa-se esquematicamente a 
constituição de um digestor deste tipo. 
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FIG. 24 


Digestor com fases cineticamente separadas, segundo 
Gosh, S. et al. [36] 


4.3 — Filtros anaeróbios 


Quando as matérias-primas não apresentam subs- 
tências em suspensão (caso típico dos efluentes indus- 
triais) têm-se utilizado com bons resultados filtros 
anaeróbios (vejam-se Scammell, W. [8], Chian, E. 5. k. 
[37] e Dewalle, F. B. et al. [38]). Nestes casos a 
matéria-prima tem de atravessar os filtros (constituídos 
por fragmentos de quartzo, lâminas plásticas,...) onde 
está provavelmente fixada a biomassa activa. 


Na figura 25 apresenta-se o caso de um filtro 
anaeróbio homogéneo, homogeneidade essa garantida 
por uma recirculação continua do efluente do filtro. 


Sonsa de 
Rectcucação 


ENTRADA PA 
MATÉRIA PARA 


soLIDdos 
PECAUTADOS 


FIG. 25 


Filtro anaeróbio homogéneo, segundo Chian, 
E. S. K. et al. [37] 


4.4 — Digestores do tipo de escoamento vertical 


(up-flow]) 


No tratamento de efluentes industriais de grande 
débito, a utilização de filtros anaeróbios torna-se dificil 
dado que são facilmente colmatados. Assim surgiram 
recentemente (ver Heertjes, P. M. et al. [39]) diges- 
tores do tipo de escoamento vertical [up-flow), em 
que o efluente a tratar é introduzido no fundo de um 
reservatório aberto de grandes dimensões. Dentro deste 
reservatório um sistema de captura de gases permite 
recolher o biogás. É então facilitada a decantação das 
partículas sólidas (figura 26) e a saída do efluente 
tratado faz-se pela zona superior do reservatório. 
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AR ATnOSEÉRICO 
EFLUEMTE 


biogas 


ALIMENTAÇÃO 


FIG. 26 


Esquema tipo de um digestor de escoamento vertical, segundo Horton, R. [33] 


4.5 — Digestores do tipo de escoamento sólido 
(plug-flow ) 


No tratamento de resíduos das explorações agri- 
colas-animais, a matéria-prima apresenta-se sob a for- 
ma pastosa, o que dificulta a sua homogeneização. 
Assim, Fry, L. J. [40] construiu um digestor (de ali- 
mentação intermitente) do tipo de escoamento sólido 
fplug-flow), em que a matéria a tratar é lentamente 
arrastada ao longo do digestor (figura 27). 


EFLUENTE 


SECÇÃO AM 


CABO DE ARRASTE DA RÊDE. 


5.1 — O biogás: características e segurança no arma- 
zenamento 


O biogás é uma mistura de gases na qual cerca 
de 60 % é metano (CH,). Os restantes 40 % são cons- 
tituídos essencialmente por anidrido carbónico (CO,) 
e outros gases, tais como, o gás sulfídrico (H,S) e o 
hidrogénio (H,), ainda que estes últimos em muito 


pequena percentagem. 


DEPOSITO PE 
HISTURA 


REDE PARA QUEBRA 
DA ESPUMA 

- E Er ERÊDE 
CILINDRO DE FLUTUAÇÃO 


FIG. 27 


Digestor do tipo de deslocamento sólido (plug-flow), segundo Horton, R, [33] 


5 -— OS PRODUTOS DA DIGESTÃO ANAERÓBIA 


Como se indicou atrás (figura 2) a digestão anaeró- 
bia pode considerar-se constituída por duas fases: uma 
ácida e outra metanogénica. Segundo Cohen, À. et al. 
[23] os produtos fundamentais da fase ácida são o 
anidrido carbónico, o hidrogénio e os ácidos voláteis 
(acético, propiónico, butírico, fórmico,...), enquanto 
que os da fase metanogénica são o metano e o anidrido 
carbónico obtidos a partir dos ácidos e do hidrogénio 
formados durante a fase ácida. 

Vejamos as características fundamentais e proble- 
mas de armazenamento do biogás bem como do resi- 
duo sólido-líquido utilizado como fertilizante. 
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Nos quadros Il e Ill, resumem-se segundo Whea- 
tley, B. |. [24] as propriedades químicas, físicas e 
fisiológicas mais importantes destes gases. 

Dohne, E. [25] apresenta alguns dos problemas 


relativos ao armazenamento do biogás. Entre eles sa- 
lientemos: 


— Perigo de explosão e incêndio — a mistura 
metano-ar só é explosiva na gama de concentra- 
ções: 5-15 %. No gráfico junto (figura 28), segun- 
do Mosey, F. E. [26], apresenta-se o domínio de 
explosão da mistura metano-ar. 
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OUDARO Il 
Proprisdades químicas e físicas do biogás produzido na digestão anaeróbia 


Biogás tipico 


Propriedade (60 94, CH, — 40% CO;) 


| % volume 100 


Valor energético 
| (Kcal/]) 


Gama de explosão 
(% vol. em ar) 


Densidade PTN 
(9/1) 


Peso específico 
(relativo ar) 


Temperatura 
crítica (ºC) 


Pressão crítica 
(bar) 


Odor 


nenhum 


E» 


ovos 
podres 


ogas 


Ss 
(a) 
E 
[a] 
E 
8 
fu 
Ê 
3 
$ 
ea 
FIG. 28 
Á volume de biogás no ar Mistura explosiva biogás-ar [26] 


QUADRO III 
Propriedades fisiológicas do biogás produzido na digestão anaeróbia 


Periodo 
de exposição 


0) 


Conc. no ar atmosférico 


Gás componente do biogás Efeitos Ffislológicos 


(o volume) 


ASFIXIANTE 
500 000 | dores de cabeça, não tóxico 


ASFIXIANTE 


inofensivo 
resp. dificil 
sonolência, dores de cabeça 
asfixia 

pode ser fatal 


VENENOSO 

irritação dos olhos e do nariz 
dores de cabeça, vertigens 
náuseas, excitação, insóneas 
inconsciência, morte 


(7 — Considera-se a exposição de um adulto humano de peso médio 68 Kg, ou outro animal de peso equivalente, Animais de peso infe- 
rior são afectados mais rapidamente, e animais mais pesados são afectados mais lentamente suportando pois niveis de concentração 
superiores. 
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— Toxicidade — o principal problema é o da 
existência de gás sulfídrico (só perceptível em 
baixas concentrações), que sendo mais denso que 
o biogás ou o ar atmosférico, tende a acumular-se 
nas zonas mais baixas. É portanto conveniente que 
as instalações de armazenagem sejam bem areja- 
das. O anidrido carbónico ainda que asfixiante não 
é tóxico e por isso não lavanta problemas. 

— Corrosão — o gás sulfídrico é altamente 
corrosivo mas só em presença do ar, logo não hã 
perigo de corrosão no interior das instalações de 
digestão. Deve apenas, tanto quanto possível, evi- 
tar-se o contacto do biogás com metais tais como 
o cobre ou o latão. 


Quanto à armazenagem, Dohne, E. [25] recomen- 
da o pré-tratamento do biogás com hidróxido férrico 
para eliminação do gás sulfídrico (ver figura 29). Aten- 
dendo a que a liquefacção é dificil, considera para o 
armazenamento: 


— a compressão a altas pressões (200-300 bar); 

-—— a compressão a médias pressões (10-20 bar); 

— o depósito a baixas pressões (do gás seco ou 
húmido) em gasómetros (a pressões da ordem 
dos 50 mbar); e, 

-—— a solubilidade em hidrocarbonetos líquidos (no- 
meadamente no propano). 


A utilização de altas pressões para a compressão 
do gás, ainda que viável exige compressores especiais 
e a energia consumida é da ordem de grandeza da 
armazenada, daí o não ter grande interesse a utilização 
deste método. 

Mesmo a utilização da compressão a médias pres- 
sões não é corrente dada a necessidade de utilizar 
compressores construídos com materiais anticorrosivos 
de custo relativamente elevado. 


biogás isento 
4 de H,5 


retentores 


válvulas de . 
balde plástico , pega e a 
4 Com O Fundo furado 
biogás do 
digestor 
FIG. 29 
Dispositivo de absorpção do H,;S, segundo Dohne, E, [25] 
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FIG. 30 


Armazenamento do biogás a média pressão [25] 


O armazenamento em gasómetros (figura 31) a 
baixa pressão, é o processo correntemente utilizado. 
Actualmente começam contudo a utilizar-se balões de 
gás comprimidos por cargas (ver figuras 32, 33 e 34), 
e mesmo a construirem-se digestores cuja cobertura de 
isolamento com o ar atmosférico é uma capa plástica 
funcionando como reservatório de gás (ver figura 35). 


FIG. 31 


seesssssessit” 
; 


GAS 


RA cn /| 


FIG. 32 


Utilização de balões para gás [25], Em construções do 
tipo faça você mesmo utilizam-se actualmente câmaras 
de ar de camião ligadas entre si 


FIG. 33 
Exemplo de utilização de um balão com uma 
câmara única [25] 
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FIG. 34 


Utilização de uma membrana elástica [25] para armaze- 
namento do biogás 


nENBRANA PASTICA 


reservatório de 

matéria prima 
E Miscalhaçê 

alimento 3) 


FIG. 35 


Digestor com cobertura deformável, funcionando como reservatório 
para gás [25] 


A solubilidade em hidrocarbonetos líquidos, sendo 
um processo que exige uma técnica altamente especia- 
lizada e instalações de custo elevado não tem sido 
utilizada. 


5 — O FERTILIZANTE 


O resíduo sólido-líquido da digestão anaeróbia 
pode ser utilizado como fertilizante (adubo não-sintéti- 
co). Efectivamente segundo a SEDE [27], os estrumes 
a utilizar como adubos deverão ter um índice C/N de 
15/1 ou inferior. Ora os estrumes utilizados na diges- 
tão anaeróbia possuem em geral valores da ordem de 
25 a 30/1, e o tratamento anaeróbio permite reduzir 
estes valores aos aconselhados para os adubos. 

Segundo Buhr, H. O. et al. [28] o manual da 
WPCF (water pollution control federation) recomenda 
que, para a destruição dos organismos patogénicos 
eventualmente presentes, estes resíduos sejam tratados: 


-— termicamente durante uma hora a 5B7ºC; ou, 


—— secos à temperatura ambiente durante 12 a 15 
meses. 
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Dal a vantagem da utilização de uma biomassa 
termófila dado que a residência (10 dias a 50ºC) das 
matérias no digestor é suficiente para a eliminação 
desses organismos patogénicos. 

Todavia outros problemas poderão ser postos no 
que se refere à utilização destes resíduos como fTferti- 
lizantes, nomeadamente: 


—= modificação do pH dos solos; e, 
— modificação da textura dos terrenos. 


Nestes domínios poucos estudos têm sido efec- 
tuados, não se tendo encontrado referências conclu- 
dentes. 


& — CONCLUSÕES 


Não se pretende com este artigo analisar exaus- 
tivamente toda uma vasta tecnologia em rápida expan- 
são, mas apenas abordar descritivamente a arquitec- 
tura dos digestores e o processo da digestão anaeróbia. 
Pretendeu-se fornecer aos interessados neste domínio 
(metanogénese) uma bibliografia suficiente para iniciar 
os seus estudos, divulgando e experimentando se possi- 
vel a construção de protótipos. 

Realizou-se pois uma análise descritiva da tecno- 
logia e processo da metanogénese. Analisamos em 
particular os aspectos quantitativo e qualitativo da ca- 
recterização das matérias-primas, a noção de biomassa 
e a arquitectura dos reactores biológicos. Fez-se tam- 
bém uma análise dos produtos da digestão anaeróbia, 
com particular incidência sobre as características e 
segurança no armazenamento do biogás. 
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VOLUME XLI 


Anaerobic digestion control using the phase-plane analysis 
(With application to methanogenesis) 


RESUMO 


Estuda-se neste artigo a utilização da análise no 
plano de fases no controlo da actividade biológica de 
digestores homogéneos, sem que se faça a lavagem da 
biomassa presente, ou se entre em «estado de choque». 


O — NOMENCLATURE (') 

A —biological activity (mg biomass . 7. dd?) 

A. — biological generation activity (mg biomass . 
rn 

A, — potencial biological generation activity of in- 
coming nutrient (mg biomass . |!) 

b  — specific death rate (d”) 

f — relative inflow of influent 


fi(L) — inhibition function for specific growth rate 
model (mg inhibitory . 1) 

F  —inflow of influent (and/or of effluent) in the 
digester (1. d7) 


k, | — inhibition constant for biomass, expressed as 
potential biomass production (mg inhibitory.. 
mg biomass . mg substrate”! . |!) 

ks — saturation constant for biomass, expressed 
as potential biomass production (mg bio- 
mass . 1) 

kj — retention coefficient for the substrate in the 
reactor (k,=0, total retention) 

ks —fixation coefficient for the biomass in the 
reactor (k,=0, total fixation) 

K, — inhibition constant for biomass (mg inhibi- 
tory . 1-7) 

Ks — saturation constant for biomass (mg subs- 
trate . 177) 

| — gubstrate concentration, expressed as poten- 
tial biomass production (mg biomass . [?) 

L  — substrate concentration (mg substrate . 7) 

R$ — net transformation rate of B into € by bio- 


logical activity (mg C . mg B'?) 


SILVESTRE D. ANTUNES (*) 
MICHEL J. INSTALLE (*) 


ABSTRACT 


Phase-plane analysis of homogeneous digester sys- 
tems is used to give information on how to control the 
biological activity without washing out the reactor or 
reaching a «shock state», 


S  — active biomass concentration (mg biomass . 
It 

t — time (d) 

Vp — reactor liquid volume (1) 

Yg — yield coefficient for the transformation of B 
into C by biological activity (mg C.mg 
BI. d) 

e  — relative growth rate for blomass 

E  — specific growth rate for biomass (d”?) 

U -—— maximum specific growth rate for biomass 
(d") 

Tt  — relative time, generation, 


1 — INTRODUCTION 


The anaerobic digestion process has been modelled 
in 1972 by Graef, S. P. [1] considering one bacterial po- 
pulation in the biological phase, Later on, in 1977 Hill, 
D. T. and Barth, C. L. [2] considered two types of popu- 
lations with differenciated functions, Actually Bryant, M. P. 
[3], Pfeffer, J, T. [4] and Zeikus, J, T, [5] are considering 
the existence of others types of populations, for which 
parameters are not yet available, 

In this paper, we investigate on the evolution of 
homogeneous population in the phase-plane, This research 
is applied to the study of the dynamics of the acid-formers 
and the methane-formers considered by Hill, D. T, and 
Barth, C. L. [2]. 


(*) Université Catholique de Louvain, — Lab. d'Automatique et Analyse de Systêmes, 
Bátiment Maxwell — Place du Levant — 1348 — Louvain-la-Neuve — Bélgique. 
() For simplification of the notations, the argument time will be omitted for variables depending on time. 


Manuscrito recebido para publicação em 4/1/B0. 
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2 — HYPOTHESES 


Let us consider an anaerobic digester (biological 
reactor-figure 1) for which the characteristics of influent, 
effluent and reactor are assumed to be: 


— INFLUENT: 


the influent in liquid phase (homogeneous) is 
introduced in the reactor at a rate F, with concen- 
trations: 


L, in substrate, and 


S, in active population (we will consider 
later S,=0). 


035 phase 


influent : F,sS effluent : EF, S.. Lj 
dinester : LE ky a ka 
Viquid phase : 5,L 
biomass : Lo, K [o 
Ei Pe id Te Fa Gt) 
FIGURE 1 
Anaerobic digester 
— DIGESTER: 


the liquid phase in the digester has a volume Vp, 
and its mean (') concentrations in substrate and active 
population are respectively: L and 8, 


The digester geometry and heterogeneity are charac- 
terised by two coefficients k; and k;, retention coefficient 
for the substrate and fixation coefficient for the active 
population. 


The active biological phase is characterised by: 


— the specific growth rate W, defined as in the 
Haldane-Monod model [6], augmented with one inhi- 
bition function f(L): 


da 


EL 
= RILFEDK aii 


— the specific death rate b, considered as cons- 
tant; 


— the biological activity, defined as: 


A=u4.sS (2.2) 


() Spatial average. 
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— the yield coefficients for the transformation of 


Cc 
one product B into C by biological activity: Yp, So 
that the net transformation rate of B into C, is defi- 
ned by: 


p= VA (2.3) 
— EFFLUENT: 


the effluent in liquid phase (homogeneous) leaves 
the reactor at the rate F, and has concentrations: 


in substrate: Ly=k:.L (2.4) 
in active biomass: Si = ko . 8 (2.5) 


We neglect the gas phase in the model as it may 
be considered as a function of the biological activity, 
and it does not introduce an important change in the 
volume of the liquid phase. Graef, S. P. [1] Hill, D. T. 
and Barth, C. L. [2] and others considered that the yield 
coefficients for gas production from biological activity 
are constants, So biological activity is enough to define 
gas production. 


3 — PROCESS DYNAMICS 
The process dynamics is described by the balance 
equations for substrate utilisation and population dynamics 


for biomass: 


substrate utilisation: 


dL o F (= kl) 1 Ss 
dt Vp ne va xá 
population dynamics: (3.1) 
dS F 
—— = (S,—ko.S)+(4—b) S 
E” Vo ( 2.8) + (1 ) 


Let us consider: 


He 
Es ne relative growth rate for biomass (3.2) 


| =Y?.L = substrate concentration, expressed as po- 


tential biomass production (3.3) 

k=Y2.K, (3.4) 

k=Yç5.K; (3.5) 

Tt=b.t= relative time (3.6) 

É = = relative inflow of influent (3.7) 
D 
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Then 3.1 may be written: 


dl 


mt em FO, | —k «=E P-— 
dr (Io 1 ) 
(3.8) 
ds 
— = ft. S,+(e—1l—k,.f). 5 


And we may define biological generation activity as: 
A =E.8 (3.9) 


The incoming nutrient introduces a potential biologi- 
cal generation activity, defined by: 


A =(h—k.D.f (3.10) 


we can neglect (if the digester is working well) the term 
for the effluent, so in a first approach we can say that 
the nutrient inflow introduces a potential biological gene- 
ration activity: 


Angsh t (311) 
Às usual we do not introduce biomass with the In- 
fluent (S, = 0) and we can write 3.8 as: 


dl 


— = f(lh—-k; .)- E. 5 
F (lo j = 1 


(3.12) 
ds 


.=(—j—k., .5 
= ( 2.8) 


4 — PHASE-PLANE ANALYSIS 
4.1 — Phase-Plane: |-S 


The state of the process may be described by the 
state variables: concentration in substrate (|) and biomass 
concentration (S) in the digester, 

Hence we may represent the evolution of the state 
of the system in the plane |-S (phase-plane, see figure 2). 

Two cases are considered in our work depending on 
the type of inhibiting function F(L): 


1 — no inhibition by the substrate: 


E(L) = c (constant) (4.1.1) 
with c independant of L (see figure 3), 
ond — inhibition by the substrate: 
f(L) = L (4.1.2) 
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S S+aS B 
FIGURE 2 


Phase-plane evolution curve 


In this case, the substrate or a linear function of It 
(for example the non-ionized substrate at constant pH,...) 
acts as an inhibiting agent for the biomass growth (see 
figure 3). 


inhibiting substrate 


non-inhibition by the substrate 


FIGURE 3 


Phase-plane supplied with relative growth coefficient for biomass 


Supplying (as shown in figure 3) the phase-plane with 
the relative growth coefficient for biomass will give us 
some suplementary information on population dynamics 
for free systems as we shall see in 8 4.3. 


4.2 — Equal activity curves 


We may consider equal activity curves and represent 
them in the phase-plane: I-S. 
From 3.9, 3.2 and 2.1 we have: 


c.s 
Á DO TT L 
É TAKAI +E(LIK, presto) 


IF the inhibition function is known, it ia possible to 
compute I(S) with A, kept constant, 
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4.3 — Free system dynamics 


Let us consider a batch system, that is a not fed 
system. Hence: 


f=0 (4.3.1) 
then the system dynamics 3.12 becomes: 
EM = — E . S 
dt 

(4.3.2) 

ds 
mm == [E me 1), 5 
dt 


and we can see that in the phase-plane the trajectory 
for the free system evolution is such that: 


dl E 
am 4.3.3) 
ds 1—E 


It is easy to show that 4.3.3 consists in a family of 
parallel curves in the phase-plane, generated by horizon- 
tal translation (see figures 4 and 5). 


FIGURE 4 


Free system dynamies for a population with no inhibition by the 
substrate 


Note that the evolution of the biomass concentration 
depends on the relative growth coefficient: 


— for E>1 
with time, 


the biomass concentration increases 


— for E<] 1 the biomass concentration decreases 
with time. 


The maximum of the generation activity occurs also 
at the maximum of the relative growth coefficient for a 
certain biomass concentration, so substrate concentrations 
greater than those giving maximum relative growth coeffi- 
cient should be avoided, as they result in lower activity 
systems (by inhibition), see figures 6 and 7, 
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FIGURE 5 


Free system dynamics for a population with Inhibition by the 
substrate 


FIGURE 6 


Equal generation activity curves for a population with no inhibitian 
by the substrate (there are no limit to substrate concentration) 


f NANA NS, a od, " 
CRM e RR Aga? Ag 

NO , j Es WE 

NORA (REGION TO BE AVDIDED ) Nas 

, - | A42 493 
| [o - | 
» » 


Ags> Au 


FIGURE 7 


Equal generation activity curves for a population with inhibition by 
substrate 


44 — Forced system dynamics 


Constant and periodic feed-rates have been studied 
and we find that if there exists active biomass in the 
reactor, the state of the reactor will tend towards a sta- 
tionary value, which depends on two limiting conditions 
for influent: 


— one feed-rate condition — limiting the value of the 
feed-rate to a maximum, depending on the inhibiting func- 
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ton (QUANTITY CONDITION FOR THE EFLUENT). So, 
if there are no substrate inhibition (inhibiting function of 
type 4.1.1), we have: 


. | 
| RS RR k (4.4.1 
io K+YE. e : a) 


and, if there are inhibition by the substrate (inhibiting 
function of type 4.1.2), we have: 


(e fo da 
1 +2 vkk 


— one substrate condition — limiting the substrate 
concentration of the influent to a minimum value (QUALITY 
CONDITION FOR THE INFLUENT) depending on the 
inhibiting function. So, if there are no substrate inhibition, 
we have: 


ko. k (1+Hk. 
Sad! END (4.4.3) 
Eck—Qi+ko. BD. (K+YE O) 


and if there are inhibition by the substrate: 


| PO Rr Aa 
[1 Viera) | 
k Ki NIttg of | 


and on the initial state of the reactor. 

If these two limiting conditions for the influent are 
verified, depending on the initial state of the reactor: 

— a stable node (N) for the working digester may 
be found; or, 

— an other stable node (WO) for the non-working 
digester (by washing out the biomass) may be found (see 
figure B). 


| 


FIGURE 8 


Trajectory evolution in the phase-plane of a forced system (influent 
potential biological generation activity: A, constant) for different 
initial state of the reactor 
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If one of the limiting conditions for the influent is 
not verified the biomass will be washed out, the digester 
will not work, and the stable node (WO) of washing out 
will be found. 

For a fixed initial state of the digester (S=Suyx and 
l=lyx), there exists one unstable node (UN), but any 
small perturbation will move its state to one of the stable 
nodes: the working one (N) or the non-working (WO). 

Furthermore, given a certain initial state for the 
digester, transient trajectories to a stationary state have 
been studied and table | summarises the conditions for 
the forced or free evolution, 


5 — APPLICATIONS TO METHANOGENESIS 
5.1 — Acid-formers reactor 


Using Hill, D. T. and Barth, C. L. [2] parameters for 
acid-formers, a simulation of an acid-former reactor has 
been implemented on a PDP 11/20 computer, The results 
are shown in figures 9, 10 and 11. 

Figure 9, represents the free system dynamics and 
equal activity curves for acid-formers, considering cons- 
tant inhibition, independant of substrate concentration, 


TABLE | 


Typical scheme of evolution after influent change 


Initial system state: |, S 
Initial system PENA Aç=E Ss 


Y 
New influent characteristics: |. E 


Lam 


Potential ii activity: A, I,. F 
À — — 


| a 
Evolution like a free 
system until: A, = A 


yes 


> (Ea) 


yes P 

no Constant activity evolution 
towards a stationary node, 
if quality and quantity con- 
ditions for influent are ve- 


de 
Evolution with 
increasing | and 
decreacing A,, until 
total washing out of 
the biomass. 


rified. 


de 
Quick increase of |, 
and of A, until Ag E As 


Evolution with increasing | 
and variable A, (increase 
followed by decrease) until 
total wash out of the bio- 
mass. 


Constant activity evolution 
until stationary node, if 
quality and quantity condi- 
tions for influent are veri- 
fied, 
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FIGURE 9 


Phase-plane for acid-formers reactor: free aystem dynamics (') 


FIGURE 10 


Phase-plane for acid-formers reactor; constant feed of Influent 


Figure 10 represents the forced system evolution, in 
the case of a constant feed rate of influent, and we can 
see that: 


-—there are no unstable node, and, 
— the washing out node can only be obtained If we 


start with a system without active biomass. 
Figure 11 represents the forced system evolution in 


() Note that, for low substrate concentrations, the “1” 
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the case of a periodic feed rate of influent, and we can 
see that: 

—trajectories oscilate around mean values that are 
almost similar to these obtained with a constant feed-rate 
of influent (see figure 10), and, 

-—— a stationary state are the stable node (N) may be 
find as it was the case for a constant feed-rate of the 
influent. 


scale becomes logarithmie. 
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FIGURE 11 
Phase-plane for acid-formers reactor: periodic feed of influent 


5.2 — Methane-formers reactor As in the case of acid-formers, we can see (figure 14) 
that a periodic feed-rate of the influent gives: 
— mean-time trajectories of the same type from that 


As in 8 5.1 the same types of curves have been made of constant feed of influent (figure 13), and, 
for a methane-formers reactor, and they are represented — a stationary state around the stable node (N) may 
in figures 12, 13 and 14 respectively. be find. 
2 
f | f O E ag 
é E e ã a 
fi ê ai E dá a 


ooo e 0 qo? Loo Tao 4uio 54% É-so “qo” J% mw q 
FIGURE 12 


Phase-plane for methane-formers reactor: free system dynamics 
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However, for constant feed-rate of the influent there 
exists an unstable node (UN), and a washing out con- 
dition may occur which is essentialy depending on the 


substrate concentration. 


We can see that a certain biomass concentration 
gives us the maximum of biological activity at the maxi- 
mum relative growth coefficient, and this corresponts to 
a well defined value for the substrate concentration. 


FIGURE 13 


Phase-plane for methane-formers reactor: constant feed of iniluent 


se 


o 


FIGURE 14 
Phase-plane for methane-formers reactor: periodic feed of Influent 


6 — CONCLUSIONS 


To model a digester with only two state-variables: 
the biomass and the substrate concentration in the reactor 
is not possible. However in a well working digester and 
assuming that all the variables are constant, the phase- 
-plane method may give us the necessary insight to impro- 
ve and to control the dynamics of the biological activity. 
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Furthermore, changing the feed-rate or characteristics 
o* the influent will change the stationnary state of the 
digester (table | summarises the evolution in this case). 


Quality and quantity characteristics of the influent 
(l, and F) define the stable and unstable nodes in the 
phase-plane and table || summarises the conditions for 
nodes existence. 
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TABLE |l 
Stable and unstable nodes for forced stationary regime 


Influent characteristics: 
|, — substrate concentration, 
F — influent edi 
Y 
QUANTITY CONDITION FOR INFLUENT: 


mean F < (Emax—1)/ ko 


yes 
e 
QUALITY CONDITION FOR INFLUENT 


mean |, > Im 


yes. | no 


Ta stable nodes and 
one unstable node (the two 
stable nodes are the nodes q q 

of working system N and There is only one stable 
the washed out system node (the node of the 
WO). washed out system WO). 
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Corporation financing under falling share prices 


ABSTRACT 


The present paper is an application of control techni- 
ques to the determination of the optimal financial policy, 
for a firm whose shares' market value falls with increasing 
leverage. By the application of Pontryagin's principle in 
continuous time the optimal dynamic policies are derived 
as funtions of the values of the parameters and the 
results are discussed. 


1 — INTRODUCTION 


In Caldeira-Saraiva (1979), we presented a model to 
study the influence of the fall in share prices with leve- 
rage on the choice of an optimal gearing policy, We 
assumed that both the rate of fall of share prices and 
the maximum amount a firm could borrow per unit time 
were small, and we studied the resulting model using 
control techniques, In the present paper, we wish to lift 
the restrictions imposed on the rate of fall of share prices 
and the rate of issue of bonds, We shall use the same 
model as previously, which we shall now proceed to 
present summarily, For more details, see Caldeira-Saraiva 
(1979). 


2 — THE MODEL 


The total capital of the firm is K(t) and the book 
value of outstanding debt B(t). Thus the value of equity is 


E(t) = K(t) — B(t) 
the gearing ratio being 


B(t) 
E(t) 


g(t) = 


If the gross rate of profit is R, the rate of interest r, 
and the tax rate (1 —0) (all constants), the rate of 
accumulation of profits is 


OIRK-r B] 


Professor do Instituto Superior Técnico. 
Manuscrito recebido para publicação em 16/7/B0. 
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by F. CALDEIRA SARAIVA (*) 
RESUMO 


O presente artigo é uma aplicação de técnicas de 
controlo à determinação da poltica óptima de financia- 
mento para firmas cujas acções tem valores de mercado 
decrescentes com a taxa de endividamento. Por aplica- 
ção do Principio de Pontryagin em tempo continuo as 
politicas dinâmicas óptimas são obtidas como funções 
dos valores dos parâmetros, sendo os resultados finais 
discutidos. 


o* which a percentage 
Y € [01] (1) 


wil! be paid out as dividends. If the firm issues shares 
through underwriters, we shall suppose the rate of issue 
is at maximum a given multiple of the going rate of pro- 
fits, i.e. 


rate of issue = n [R K—r B] 
with 
0O<n<N (2) 


and that a percentage, proportional to gearing, is charged 
by the underwriters to cover their risks, 
The proceedings of the issue for the firm are then 


(1-2 9) n [R K—r B] 


and the rate of increase of the equity will be 


dE 
a = 1 0+(1-ag)n] [RE+B)-—rB]) (3) 


The firm can also obtain finance by borrowing. The 
amount of new bonds a firm can issue per unit time will, 
similarly to shares, be at most a certain number of times 
U, the going rate of profits, i.e. 


dB 
= =b IRE+B)—rB (4) 
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with b <. U. The firm can also repurchase bonds out of 
its retained earnings and the proceedings of the sale of 
new shares, so 

U>b>-[1—-9ê+(1—a gn] (5) 


The objective of the firm is the maximization of the 
owner's discounted receipts: 


T 
Vío) = | <(d(t) — dt) e tt dt + [E(T) — EL (M)] er 


“o O 


where d(t) is the dividend at time t, T the (supposedly 
constant) income tax rate for shareholders and d.(t) and 
E (t) represent respectively the dividends and equity of 
new owners (owners different from the owners at time t). 
Adopting the procedure introduced by Krouse (1972) and 
followed e.g. in Krouse (1973), Krouse and Lee (1973), 
Senchak (1975), Perrakis (1976) etc... we can equate the 
discounted value of shares placed with underwriters with 
the discounted value of dividends and final equity, The 
objective then becomes: 


cm. 


V(o) = | (ty 0—n) [R(E+B)—r B] et dt+E(T) e ET 
“O 
(6) 


which we want to maximize, subject to the restrictions on 
the control variables (1), (2) and (5), the development 
laws of the state variables (3) and (4) and the restriction 


B>0 (7) 


3 — SYNTHESIS OF OPTIMAL SOLUTIONS 


Taking adjoint variables LVp=9p e *t and Y=9a CH, 
and a shadow variable Ag et! for restriction (7) the 
Hamiltonian of the system is 


Jt= | 70—n)] IRE+B)—r B] et + 
+ (gp e*) [(1—7) 0 + (1—a 9) n] [R(E+B)—r B]+ 
+ (Pp e) b [R(E+B)—r Bl+(Ap et) B 


Le] 
Putting H = df e*t and representing by X the deri- 


vative of variable X with respect to C=T-—t (backwards 
time), the equations of the system become 


E) da B 
(E) = ([1—90+(1—09) n] g—b) 


do E 
[R+(R—r) 9] (8) 
o | a St â 
pe +k Qp= e! e ah = akt ari H 
q dt dE dE 


=(7780-n) R+oell(l—1)8+(1—a 9) n] R+ 
+eng [R+(R—r) 9]) + 


db 
+ és [SR=IR+ (A 9] 9 Tr] (9) 
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=(T70-—n) (R—-)+9; (MI—r) 0 + 


é d H 
Pptk Po= 5 
+(I—a 9) n] (R—-)—a n [R+A(R—r gl) + 


| db 
+Pa |pe-o+R+n=o g] | +hg 


am fa É )- (xy 0—n) [RA(R—N) gl + 
+(I— 90 +(I—a 9) n) [RAH(R>—N) gl—k) w 
(1) 


where 7, n and b are chosen to maximize 


+IR (E+B)—r BJ (y (t—eg)+n [(l—a 9) qe—I]+ 
+b Pp) 


or 
H=H+IRE+B)-r BI (1 S, +n S,+b 8) (12) 
where H is independent of the control variables and 
ds 
Sn=(1—0 9) Pp—1 
Sp=Pp 


when 7=0, Pg is one and Pg is zero by the transver- 
sality condition, so 


Sy=t—l 
9n=—€ 9 
9,=0 


The income tax rate T can only be less or equal 
to 1. If v<21 then S,<0 IF t=1 then >, =0 but 


from (9) 


S,(0)= — PplO) 
=—(0R-k—a ng [R+A(R—r) 9] 
And, as by hypothesis O R > k, (or else the business 


after taxes would not be worthwhile) it follows that 
S, (0) < 0. So, in any case the optimal policy for divi- 


dends is 7=0. 
If gq, > O then S,(0) will be negative so the optimal 
n will be zero; S,(0) but S,=8 (R—r), so S(0) > 0 


and therefore, the optimal b is b=U. 
If g,=0 then S,=8S,=0 and 


| 


o O o 
Sn= "ug Pe+(i—a 9) Pe 
[a o 
=abR+(0 R—k) 


S, =0(R—r)—a n R+Ag 
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Let us first suppose that 
GR—-)-a NR>O 


Then S, will be positive, so at 7=0t we shall have 


S, > 0 and thus b=U. It follows then that 9(0)= —UR 
and thus at 9=0*, q < O thus violating constraint (7). 
So, there can be no backward trajectories starting at 
g,=0, 1.e. all real-time trajectories end with g, > 0 and 


A A A 
controls at Y=0, n=0 and b=U, Equations (8), (9) and 
(10) yield 


9 =(8 g—U) [R+(R—P) gl (13) 
Pe=(0 R—Kk) QE+U R (14) 
Pp=UR—r) Pe +HUR—N—K] 9p (15) 


and, integrating 


U[R+(R—r) go]+R (8 g,—U) elf R+U(R-n]o 


q A 
IR+HR—A) q]—(R—9 (8 go—U) e É PACTO 
(16) 
o Rel B+UGOOO uyp-r) ko " 
id ORAU (RO) o eo 
0 (R— [0 R+U(R-"] —k 
” A)  (RAUR-O | to 


“0 R+U (R-r) 


From equations (17) and (18) we see that q; and Pp 
will remain positive. Therefore, from (14) it follows that 
PE will be positive and so Pp will remain above 7 and 
there will be no switch in 7. As S,=?p, there can also 
be no switch in b. The only switch that can thus occur is 
if Si=(l—a 9) qg—l changes sign. This will happen 
for a value of q equal to 


so (1- ) (19) 
Osw = ao Pe 


The adjoint variable q, is, as we said above, a func- 
tion of q which takes the value 1 at c=0 and increases 
to co as € -—> co. Thus go) is a function taking the 
value O for € = O and increasing to 1/x as T —s 00. 

In Caldeira-Saraiva (1979) we studied thoroughly 
the case when 1/x > U/0 which is obtained for small 
values of both q and U. We shall now, as said in the 
beginning, concern ourselves with the case 1/u < U/B. 

If g, > U/8 then, we see from (13) that g will increa- 
se It will therefore be always greater than U/0, thus 
greater than 1/x and a fortiori greater than gu, for the 
same 9. Thus there will be no switches. The same will 
happen if g,=U/8. If q, < U/8 then g will decrease ten- 
ding to zero In a finite time. The descending patch will 
therefore cross the rising switching curve g.,, and there 


TÉCNICA 459 


will be a switch in control ffom n=0 to n=N. The system 
will then evolve under: 


9 = ([0+(I—a 9) N] g—U) [R+(R>0) 9] (20) 
dm = NR + ([8+(I— a 9) N]) R42Ng [R+ 
+(R—r) gl—k) q; +UR pp (21) 
Pp =— NR—o)+f[e+(I—a 9) N] (R—9)—a N[R+ 
+(R—rgl) 0; +IU(R—)—k] qp (22) 


Equation (20) can also be written as 


q = [R+(R—1) 9] [-=N a g?+(8+N) 9—U] (23) 
The second degree polynomial has roots 


A 0+N + /(0+N)—4 & UN 
=" 2aN 


Whether these roots are complex or real depends 
on whether the discriminant 


(0+N)º—4 À UN (24) 


is negative or not. This discriminant is itself a second 
degree polynomial in N of discriminant 4=16 a” 8 U 
(U/0—1/a) As we are assuming that U/0> 1/x, A is 
always positive and therefore (24) has always two real 
roots. Let them be N, and No. Writing (24) as 


U 1 
N—2 qa 6 [2 p= N4+ 02 
Ê o 


shows us that N,+N; and N;N; are positive, so N, and 
No will both be positive and we can assume, without 
loss of generality, that N, >= N5. Putting x=0/N we have: 


+ x V(1 + x)?—4a (U/0)x 
| 2a | 


912 00 = 


and, differentiating 


(3:--) 
da Es 91,2 — q 


dx a V (1 +x)2?— 4 a (U/0) x 


When x=0 (i.e. when N > 50), g1=1/a and g,=0. 
As 1/u< U/0, it follows that dg,/dx is negative and 


A A 
dgo/dx is positive. As x increases from zero, g, will 
decrease thus continuing to be less than U/0 and 


g> will increase but, being less than g;, will stay under 


U/8. 
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When x becomes 0/N, the discriminant will be zero 


E " 1+8/N, dg 
and g, and g; will be a with q e 68 and 
dg 
Te = +00. 


When x is between Q/N, and 0/N, (ie, when N is 
between the roots), the discriminant (24) is negative, 


and q, and q are complex. For x=0/N,, 9 and 9 will 


1 +0/N; 


both be equal to Sd Replacing N; by its value, 
A 


ER (Pl 


or, multiplying the numerator and denominator by 


U / U U | Tira ca ds E 
em ci im z - — | , and simplifying: 
0 V 4 


and, multiplying again both terms by the same factor 


À mB /U(U 10 


which shows that RP and g are greater than e 


Then, dg;/dx will be positive and dg/dx will be 
negative, so g; will rise and g will decrease as x increa- 


ses. Às x-> 00, 9; will also tend to vo and > will tend 


o 4 
“a 


A A 
g, and g, with x. 


In the figure below, we see the behaviour of 


ah ole 


FIG. 1 


ga 
We saw above that for g, > U/8 there would be 
no switches and n would always be zero. Switches can 
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only oceur for g < 1/x, when the descending path with 
n=0 crosses the switching line. 

Let us now go back to the behaviour of g after the 
switch, described by equations (20) or (23). IF N<N, 
(ie. x > 9/N,) the two roots of the second degree poly- 
nomial will either be complex or greater than U/0 and 


therefore, than the going g. It follows then that g will 
always be negative. Trajectories will then be as in fi- 
gure 2, 


FIG. 2 


A A 
E N>N, there will be real roots gy, 92 in the 
interval (0, 1/0); g will be positive if ge (92, 9) and 
negative if q is outside this interval. The values a and 


Fa will be asymptotes. 

The switching curve, given by (19) has a slope 
which, as we see from (17), is totally independent of N. 
Between 3 and õ, the slope of trajectories after swit- 
ching is given by (23) and depends on N. Differentiating 
we have 


3g 


NC IR+(R—N) gl (I—-ag) g>0 


À o] 
which shows that, for given g, the value of g will increase 
with N and tend to co as N tends to c0. 
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kt follows therefore that, for a small N (such as Nº 
in the figure) the slope of the paths after switching will 
always be smaller than the slope of the switching curve. 


The paths will be as in the figure below 


FIG. 4 


However, for N sufficiently large (such as Nº in 
figure 3) there will be two points g”; and g”; for which 
the paths after switching are tangent to the switching 
surface and such that for any g in between, the switching 
curve has a smaller slope than the trajectories after 
switching. Let A in figure 5 correspond to the point of 
the switching curve for which g=9g".. 


FIG. 5 


After switching at A, the path is tangent to the 
switching curve which means that the forward continua- 
tion of the path remains below the switching curve for 
a while, intersecting it again at B. There will, therefore, 
exist a dispersal surface between À and E (the point 
for which g=9"5) such that backward trajectories that 
intersect the switching surface between B and E will, 
after switching, terminate in the dispersal surface, at the 
same point that backward trajectories that would have 
switched between A and E. 


immediately after the switch, as we saw above, the 


controls will be Y=0, n=N and b=U. We must now 
investigate if there will be another change in controls. 
This will happen if there is a change in sign in 
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9, =T—P 


9n sa (1 — g) Pe—1 
or 
Sp = Pp 


Writing equation (21) as 


pe=N RI(l—a 9) pe 1]+(0 R—k+a Ng [R+ 
+(R—N gl) 9g+U R os 
we see that as long as S, and S, remain positive, Pr 
will increase and therefore S, will continue to be ne- 
A 
gative and Y to be zero. 
If q decreases, as happens in figure 2 and below g2 


and above g, in figure 4, then S,=(1—0 9) pp—1 will 
necessarily increase and will thus remain positive. But, 
if gq is between g; and gn g will increase and we must 


differentiate Sn to find out whether there will be a 
switch in N. 


d 
S= (1 — ni 
n do [(1 (x g) Pr I ] 


Le] [k] 
=-9gY+HIi—a 9) PE 


Replacing É and Pr by the values given by equa- 
tions (20) and (21) and simplifying we obtain: 


S,=(0(U—0 9) IRH(R—9 gl+(1—0 9) (O RH) qr + 
+NR [(l—a 9) pp; —1+NR (1—a 9) Pg 


Since q will be below g: and g: is lower than 1/% 
and, therefore than U/0, it follows that aslong as QPg, 
S, and Pp remain positive, S, will increase and will 
therefore continue to be positive. So n cannot be the 
first control to switch. 

Finally 


S,=05=N(R—r) [(1— 9) Ep HUR=N-K] Pp + 
+(0 (R—-9)—a NIR+(R—N gl) qe 


We are still assuming that (R—rn)—a NR >0. If 
(R—)—a NIR+(R—r) 9] is also positive for all 
g<1/x Pp cannot change sign, which means that 
there will be no more switches. The final picture is then 
as shown in figures 2, 4 or 5. 

For a switch in b to be possible is then necessary 
that we have 


AR—rn)—a NIR+(R-n) g] <0 
” MR—r)—a NR 
2 0" aNAA 
=| 
On the other hand, in (22) the coefficient of Pp Is 


[Bri—a 9) N]) (R—)—a NIR+H(R—r) 9] 
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which is zero for 


ap (R—n)—a NRAN(R—r) 
io 2 à N(R=t) 


Todas 1 
a (0i+5) 


lf Pg becomes zero for g > q> then Pp < 0 at that 
time, so there will be a straight switch from b=U to 
b=- [D+ (1—ag) N]. 

As Yy reaches zero for lower and lower values of 
g then, before we reach a there will certainly come a 
point for which when Pg=0, PB is zero, for If we ever 
reached q! then we would have Ppg=0 and Pa > O which 
is absurd. 


When Pp=96=0, the second derivative is 


qp=a NIR+(R—r) g2—N k(R—r) — 


2wuNº(R—r) [+ (1-0 9) N] glR+(R—n 9] 
[ê+(1—a 9) NI(R—)—a NIR+(R—r) gl 


A 2 « NY4A—r) [RA(R—r) 9] 
B+((1—a 9) N] (R-r)—a NIR+(R—P) gl 


Since we must necessarily have g < g> it follows 
Í Pp 
db 


values with Pg=0, Palb=U) is surely non- negative. 


that 


> 0. Now, as Pp reached zero from positive 


Palb= — [BD +(1—a q) NJ) is necessarily a lower 
value. If it is negative then either b=U (which gives 
rise to pa) >0 o b=—[8+(1—-a 9) N] which 
causes Pp(7*) to be negative, are optima! policies. Indeed 
there must be a bt el—[8 + (1—a 9) NJ]; U] such that 


Qulb=bt)=0. This value of b gives then rise to a uni- 
versal surface !. 


Along it will shall have (from the condition Pa=0) 


N(R —r) 

Pg=—"—————————————————————, (25) 

[0+(1—a g)N] (R—)—a NIR+(R—r) g] 

o dE 
and, from g = de | —— we obtain 
dg do 

0 +(1—« 9) N] (RA—)—a NIR+(R=—n) g] 

id 2 a N(R—r? 
(x NIR+(R—r) 9P—k(A—rn) (26) 


which gives us the differential equation of the universal 
surface. 

The bynomial in braces is always positive as 
q NR>0(R—r) implies « NIR+(R—r) g/>0 R(R—y=> 


!. In the terminology of Isascs (1965). 
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> K(R—r). Then the sign of o is given by the numerator 
of the first term. This is positive for gm, so 9 
increases for q < g2, and tends to zero as g 7 q>. 

From (25) we see that Pp will rise so there will be 
no change in Y. Also 


« NIR+(R—r) g]—0 (R=") 
“[94+(1—4 9) N] (R—r—a NIR+(R-9) gl 


and, as q is necessarily less than 9; and more than Ri 
it follows that Sn will be positive. Thus there will be 
no switch in n, 


As q decreases we see from (8) that it is necessary 
that bt increases. IÉ N < N, we know that the first term 
in (20) is negative which means that any positive 
however small can be obtained with bt < U. So the 
universal surface will have an asymptote at g=92. At 
any point of the universal surface we can switch to 
either b=U or b=—[B+(1—-a 9) NJ In the first 
case we shall have Ps > 0 so Pg will become positive 
thus vindicating the choice. The equations of the system 
will then be (20), (21) and (22) and g will fall. In the 
second case Pg will become negative so the choice will 
also be optimal. The equations ruling the system will 
become. 


q =[0+(1—a 9) NJ (g+1) [R+H(R-P) gl 


Pe =— NR+qe(fO+(1—a 9) Nl+o NRIR+(R—r) gl— 
—k)—os(lo+(i—a 9) N] RHR+(R—N g «a N 9) 


o =— N(R—n)+9 ([0+(1—a 9) NI (R—-—a NIR+ 
+(R—N) g)j—es([0+(1—a 9) NIR+ 
+HIR+H(R—r) 9] &N—K) 


and we see immediately that q will rise with 0. 
o 
Since Pp can also be written as 


Pe=[(1—0 9) Pp—1] NR+(0R—k+a NR [R+ 
+(R—r) gl) o;—o; (0 +(1—a 9) N] R+IR+ 
+HR—-) g]) É N q) 


we see that as long as S,=(I—a 9) Pp—1 is positive 
and Pa negative Pp will rise, so there will be no change 
in 

The only two changes that can occur will take 
place if Pg becomes positive again or S, negative, 
If qa changes sign first b will switch to U and the 
trajectory after switching will have a lower slope than 
before thus tending to intersect other trajectories with 


TÉCNICA 459 


, =—[0+(i—a g) NJ. This means that there will be 
a dispersal surface ending the backward trajectories. 

The terms in n and b in the expression of the Hamil- 
tonean are, as we know 


P=n Satb Sh 
and, if b=-—-[08+(1—-a 9) n] we have 
P=-0 S+[S4—(1—a 9) Sj] n 


It follows then that for S,>(l—a 9) S, even if 
S,.< 0 we shall still have n=N, changing to n=0 when 
Sa<(I—a 9) Sp Now S,=(l—a 9) ppt so if g 
increases, the change must necessarily occur before 


1 
g = pre as there we would have 
Spa=—1I<(i—4 9) S,=0 


So if pg did not change sign first we shall, before 
we reach g = = change n to n=0. If q; changes sign 


after that the switched paths will have as controls b= U 


1 
and n=0. If g, though greater than = is still less than 


U 
rm g will tend to fall with g (see equation (13)) and, 


ae before we must have a dispersal surface. But if 


U 
g 2 + 9 Will rise with & and we can also have a 


straight switch. The full picture is then as shown in 
figure 6. 


FIG. 6 


We see a universal surface H | from where emerge 
downward-sloping curves with n=N, b=U and upward- 
-sloping ones with n=N, b=-— [8+(1—a 9) NJ]. These 
last ones either end in the dispersal surface E L or 
switch n to zero at switching surface L M and end in 
dispersal surface L P or by switching b to U at the 
switching surface P. R. 


E N>N, then for n=N and b=U g will be posi- 


A 4 
tive if g is between g; and g,. As we saw above the 
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universal surface lies between o and g> As we saw 
above the universal surface lies between 91 and q and 
has an asymptote at g which means that as g rises to 


q the slope necessary to keep on the universal surface 
will tend to zero. If g; falls outside the interval [92.91] 
then the universal surface can approach its asymptote 
with bt in the interval [—[8 +(1—« 9) N], U] as we 
see in figures 7, 8 or 9. 


FIG. 8 


FIG. 9 
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But if gs falls inside the interval then the universal 
surface's asymptote could only be approached for values 
of b greater than U. That means that at some point | 


where g is so near q> that b =, the universal surface 
must end and the only possible optimal trajectory is the 
one with b= — [8 +(1—a 9) N]'. The universal surface 
will be followed by a switching surface where the tra- 
jectories with b= U that originated in the universal sur- 
face will switch their control to b= — [8 +(1—a 9) N] 
as is shown in figure 10, 


FIG. 10 


In particular, if at point H the value of b necessary 
to trace the universa] surface is already too large, the 
universal surface will vanish altogether to be replaced 
by a continuation of the switching surface, as shown in 
figure 11. 


b=U 
i-299N) 
Esas | % 
ne W 
bz e | 
= | sad 
Ds a | | 
DU) 
FIG. 4 


| Having thus completed the study of the case 
Palb= —-[8+(1—-a 9) NJ]<O both for N<N, and 
N>N, let us now see what would happen if at 
the point where Pa = Pp =0 we had Pplb = — [0 + 
+(i—a q) NJ] > 0. Then, the value of b which would 
be needed at H to preserve the singular condition 
Pa = 9 = O, would be smaller than —[8 +(1—a 9) NJ] 
corresponding thus to a slope steeper than that obtained 
with b=— [8 +(1—a 9) NJ. Since this is the smallest 


1 See Kelley et al. (1967), section 3.17. 
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value of b we can select it follows that the only possible 
backward trajectory at H is the one obtained with b=U, 


As Pa is positive for b= = [8+(1—a gl] NJ at H, 
Pa will also be positive for paths that switched b from U 
to —[0+(1—a 9) N] near H; this means that q; will 
only temporarily be negative and will tend to become 
positive again. 

On the other hand as we know, the slope of the 
universal surface is smaller as it approaches its asympto- 


te so, when Pp reaches zero again, if Y,=0 we may 
have a universal surface that can be traced with admissi- 
ble values of b*. 


FIG. 12 


In figure 12 we see a case where N < N, Paths 
switch b from b=U tob=—[8+(1—a 9) NJ at the 
switching surface botween E H or terminate at dispersal 
surface E P. At | the universal surface is tangent to the 
trajectory with b=— [8 +(1—a 9) N] and it proceeds 
from there with admissible values of bt, to approach 
from below its asymptote at g = gp. 


É N>N, and g; is between g; and en then the 
universal surface will eventually come to an end when b* 
tends to go above U, as in figure 13. 


FIG. 13 
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We have finally to consider the case 
“RN > MR—r) 


As we saw above, If the opposite happens there 
AR—)—a« RN 


b itch in b f Es 
can be no switch in or g < 91 ZNR=D 


Now q, is negative and b=—[0+(1—a q) NJ] is an 
optimal policy even for some point of the line g=0. 
The analysis follows the same lines as in Caldeira- 
-Saraiva (1979) so we shall not repeat it here. In fact 
the difference between the model considered there and 


1 U 
the one we have now is the sign of — — ida But since 
x 


most of the analysis for «RN > WR-—r) indicates that 
the optimal b is negative, the value of U (and therefore 
the fact that U is greater or smaller than 8/x) is irrele- 
vant. The pictures that follow summarize thus the results: 


FIG. 14 


In Fig. 14 «RN > (0+N) (R—r) and the line g=0 
is a semi-universal surface from where trajectories with 
n=N, b=-— [0+(1—a 9) N] emerge. The trajectories, 
similarhy to what we saw above switch n to zero before 
g reaches 1/x and eventually b to U, or alternatively 
terminate at a dispersal surface. 

In figure 15, «RN < (0+N) (R—r) and the line 
g=0 is only a semi-universal surface till 


1 Nla NRZ—k(R—n] 
E E DE a io o E 
(0+N)R—k (0R—k) [(06+N) (A=9)—« NR] 


At this point == and 


S, = N[NR2—k(R es = 
ps aNRBR:—kiR— PRP RARE RAT 
ja D+ 3an R=n=aNR * 


a 


Sp increase with b and is certainly positive for 
b=U. If it is negative for b= — (08+N) then there 
will be an intermediate value of b which will start a uni- 
versal surface with the same defining equations as before 
(see equations (25) and (26) ). 


IF q lies outside the interval [92.91] (or If Ja 9: 
are complex) the universal surface will rise to its 


asymptote at q, as in figures 15 or 16. 
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FIG. 16 


If however g falis inside the interval [92.91] then, 
as before, the universal surface will stop when bt reaches 
JU and change to a switching surface as in figure 17. 


nro | pe 


te=vu 


AE 
AM NE 


4 — CONCLUSIONS 


Let us now interpret the results obtained, 
The condition 


U 1 
E na 
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means essentially that both U and « are relatively large, 
ie. that there is a good availability of credit but that 
borrowing has a veery strong depressing effect on the 
market value of shares. If N is not too large, borrowing 
at the largest rate is the policy one would expect, com- 
plemented by issuing shares If g is sufficiently low to 
make the issuing premium small and/or the time left to 
recip benefits sufficiently large to compensate for the 
premium paid to underwriters. That is indeed the situa- 
tion we observe in figures 2, 4 and 5. 

When the gearing ratio is low, B is small compared 
tô E. Then if one is borrowing at maximum rate and at 
the same time issuing shares, the effect of the borrowing 
on B is proportionally larger than the effect the share 
issuing has on E. Thus g will rise. If a is large, then for 
values of q comparatively low, the proceedings of share 
issuing are low as the percentages demanded by under- 
writers will tend to be large. If one is borrowing and 
issuing shares at maximum rates, the effect on B is 
greater and again g rises. 

However for intermediate values of q, as we see in 
figures 4 and 5, share issuing while borrowing will 
reduce the gearing ratio. In this region, it is clear from 
the considerations above that when the time left is large 
it is better to issue shares and when it is small it is 
better not. There can exist, as in figure 5, a region in 
between where it is indifferent to issue or not. This 
dispersal surface is later (i.e. for smaller values of 0) 
than the switching surface, since the issuing of shares, 
as it decreases gearing and premium associated with it, 
is itself a factor making for the atractivness of issuing 
them, thus making it possible to be undertaken nearer 
the end of the project. 


For large values of N, ie. those able to satisfy 
& NIR+(R—r) g] > MR—r) for some g < 1/a, then des- 
pite the high values of «, the large size of issues allowed 
may make share issuing atractive. So much so that as 
policy of repurchasing bonds (negative values of b) to 
decrease gearing and underwriter premiums may be indi- 
cated for large values of q and 0, as a prelude to share 
issuing. This is what we see in figures 6 to 17. 

In the last figures the size of possible issues is so 
large («NR > WR—r)) that it may even be optimal to 
reduce the gearing all the way to zero. 
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Simulador industrial de controlo de nível 


RESUMO 


Para fins didáticos construiu-se em 1976/77, no 
Centro de Automática da Universidade Técnica de Lis- 
boa (C. A. U. T. L.), um sistema para estudo do con- 
trolo de nível, utilizando sempre que possível compo- 
nentes industriais. Na nota técnica anexa descreve-se 
sumariamento este sistoma. 


1 — INTRODUÇÃO 


Talvez devido a um certo atraso no ensino do 
controlo automático aplicado na indústria, não existem 
vulgarizados na nossa Universidade modelos mais ou 
menos reais de aplicação deste, ainda que na indústria 
existam já alguns modelos normalmente utilizados na 
formação de técnicos especializados. 

O modelo construído no C. A. U.T.L. (Fig. 1 e 2) 
pretende de uma forma modesta colmatar essa falha, 
em particular em relação ao ensino superior de enge- 
nharia na Universidade de Lisboa. Nesta fase constru- 
tiva (muito se pode fazer no sentido da sua ampliação) 
tem-se a possibilidade de estudar o controlo automático 


FIG. 1 


]. R. CALDAS PINTO (*) 


ABSTRACT 


For didactical purposes, a level control system has 
been assembled at C. A. U.T. L., in 1976/77. Industrial 
components were preferably used. This technical paper 
presents a brief description of the system. 


do nível da um depósito quer actuando no caudal de 
entrada quer no de saída deste. Para o efeito está 
instalado um regulador de acção PID (Proporcional, 
Integral, Derivativo). Na Fig. 3 e 4 apresentam-se os 
dois diagramas de blocos actualmente possíveis. 

Poder-se-ão criar perturbações bem determinadas 
ao sistema através da actuação no caudal de saída do 
depósito D (Caso |) ou através da actuação no caudal 
de entrada de D, (Caso ll). 

Nas figuras 3 e 4 consideram-se os componentes 
do sistema caracterizados por um ganho constante, o 
que é localmente válido. Na determinação da função 


Circuito hidraoulico 
electrico 
" pneumatico 


FIG. 2 


(*) Assistente no |, S, T., investigador no C. A. U, T, L.— (Instituto Nacional de investigação Cientifica), na Linha de Acção 1: 


Controlo de Processos Industriais. 


Manuscrito recebido para publicação em 11/7/B0. 
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FIG. 3 


PROCESSO 


FIG. 4 


de transferência do depósito D linearizou-se a função 
Altura de líquido — Caudal de salda e daí a validade 
do modelo se restringir a pequenas variações de nível 
do depósito. 


2 — DESCRIÇÃO DO PROCESSO (v. Fig. 2) 
2.1 — Circuito da água 


O nível que se pretende controlar é o da água 
num depósito cilíndrico D. Utiliza-se um segundo de- 
pósito D, com descarregador aproximadamente linear, 
com o objectivo de introduzir um atraso no sistema. 


O ciclo da água é o seguinte: 


Entrada de água em D, regulada por uma válvula 
de actuação pneumática V., e saída para D. De D a 
água passa para um reservatório R: 


por uma válvula de actuação pneumática V., 
e/ou 


por uma válvula manual e/ou 
por uma electroválvula E V 1 


Para garantir um nível de referência constante os 
tubos de descarga do depósito D mergulham no reser- 
vatório R de nível de água constante. 

Em funcionamento normal o circuito é fechado, 
sendo a água recirculada através de V,, por uma bom- 
ba. cuja curva característica é tal que a altura mano- 
métrica de elevação pouco varia para o regime de cau- 
dais utilizado. 
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2.2 — Variáveis medidas 
2.2.1 —— Medida de nível do depósito D 


A medida de nível do depósito D faz-se por me- 
dida da pressão hidroestática no fundo daquele, utili- 
zando um manómetro com saida sob a forma de varia- 
ção de uma resistência eléctrica. 


2.2.2 — Medida da pressão da água antes da válvula 

Vo 

Para uma leitura da pressão da água à saída da 
bomba utiliza-se um manómetro estático. 


2.3 — Regulador 


Como elemento fundamental de regulação utiliza- 
-se um regulador de acções proporcional, integral e/ou 
derivada (PID), prevendo-se igualmente a regulação 
numérica (por processador ou microprocessador). 


2.4 — Adaptadores de sinal 
Em virtude do regulador ter entradas e saída em 
corrente eléctrica, utilizam-se conversores electro-pneu- 


máticos (E/P) —- um para cada válvula pneumática a 
um conversor resistência/corrente eléctrica. 


3 — VARIÁVEIS DO PROCESSO 
3.1 — Lógicas 
Detector de nivel máximo tipo boia. Utiliza-se 


como sistema de protecção e alarme. Na Fig. 5 mos- 
tra-se o circuito eléctrico que lhe corresponde. 


24V 


FIG. 5 
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3.2 — Analógicas 
3.2.1 — Variáveis de entrada 


— Comando do valor de referência do nível, Há 
possibilidade de se ter como variável de referência 
toda a função que se possa gerar electricamente. 

— Comando do caudal de saída. 

Variação deste em escalão através da electrovál- 
vula EVT ou outra função qualquer através de V., e 
respectivo conversor E/P. 

— Perturbações. 

O caudal de entrada em D, depende não só da 
abertura da válvula V., mas também da pressão a mon- 
tante e a jusante desta. Esta depende da bomba de 
elevação que a faz depender do caudal. Contudo a 
variação é pequena na zona normal de funcionamento. 


3.2.2 — Variáveis de saída 


A variável de saída é o nível de água no depósito 
D, cujo valor se pode obter por leitura da pressão no 
fundo do depósito, linearmente relacionada com a al- 
tura. 


4 — COMPONENTES 


4.1 — Regulador de acção Proporcional, 
Derivada (tipo PHILIPS PCS) 


Integral e/ou 


1 — amplificador de acção derivada 
2 — pré-amplificador 

3 — amplificador de saida 

4 — realimentação RC 

9 — medida do erro x, 

6 — medida da saída 


4.1.1 — Função de transferência 


Y 
= =K, (t+ Tess 


X 


E) 


* 
com Xp ajustável entre 2% e 400% em 12 passos 
* 
é Sd » » 2 a 30m por escalões ou 
To = 
* 
Ty » » O a 13m continuamente 


(*) Segundo o manual Philips PCS (figura 1). 
(*) Segundo o manual PHILIPS E/P. 
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sendo 


Xp = —= Xg 
* 
Ty 
T2A Tt e A=1+— 
Tx 

, w 

LE ia 


4.1.2 — Entrada: corrente eléctrica O a 20 mÃ 
Saída : » » O a 20 mÃ 


4.2 — Conversor electropneumático (tipo PHILIPS E/P) 
X 4 Y 
4.2.1 — Função de transferência 


= Sir 


4.2.2 — Característica dinâmica (*) 


HENNNRR 
DESEN dm 


= 


4.2.3 — Entrada: corrente continua 0-20 mÃ 
Saída : ar sob pressão de 0,2 a 1 bar 


4.3 ——- Conversor: resistência — corrente eléctrica (tipo 
PHILIPS) 


A Y 


4.3.1 —- Função de transferência 


Y 
x = É 


405 


4.3.2 — Entrada: 50-30-50 (0) 
Saída : O-—20 ma 


4.4 — Manômetro com saída eléctrica (tipo PHILIPS) 
x A Y 


4.4.1 — Função de transferência 


" K 
x — tur 


4.4.2 — Entrada: O a 1000 mm de H, O 
Saida : 5650-30-50 0 


4.5 — Válvulas pneumáticas e respectivos actuadores 
(tipo DRAYTON) 


4.5.1 — Equações características 


Q =7,55. 105C,Ah Q m3/s 
- x 
pa Cy mae 1 h-m de H, 0 
mo a Cy min PE Voy > Cy mar = 6 
Vp2 Pe C, max = 10 


4.5.2 — Entrada: 0,2-1 bar 
Saída : caudal 


Na Fig. 6 e 7 apresentam-se as curvas caracterís- 
ticas das válvulas utilizadas. 


4,6 — Electroválvulas (REMSA) normalmente fechadas 


5 — PROCESSO 


5.1 — Depósito com descarregador linear (D,) 


Função de transferência 


Sp — área da base do depósito 


R 
ET s+1 H — constante de proporcionalidade 


5.2 — Depósito D 


Função de transferência 


c — constante resultante da li- 
nearização de H. 

k — constante relacionada 

K com a perda de carga na 
válvula manual 

S,— área de base do depósito 


CURVA “DESLOCAMENTO DA HASTE DA VALVULA VvpI-. 


= - PRESSÃO DE ALIMENTAÇÃO. 
636 

uá 

212 

“az CO sai RR = 


00 «MOV. DESC. DA HASTE 
mw bi ASC. V1 


0,76 james nas Ê Bar 


FIG, 6 
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CURVA “DESLOCAMENTO DA HASTE DA VALVULA VP2- 
-PRESSÃO DE ALIMENTAÇÃO” 


2454 
++ MOVIMENTO DESCENDENTE DA-HASTE 


mem =» " ASCENDENTE » “ 
dr 035 Ro 0 DES. , 095 OS Bar 
FIG. 7 
CONCLUSÃO 
Este protótipo faz parte de um projecto desenvol- essencialmente a fins didáticos e nele se previu a possi- 


vido no Centro de Automática da Universidade Técnica bilidade de controlar uma segunda variável — a salini- 
de Lisboa em colaboração com o Departamento de Sis- dade da água em circulação, e também o controlo numé- 
temas (Máquinas) do |. S. T., e no âmbito da Linha de rico por computador. 

acção «Controlo de Processos Industriais». Destina-se 
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Moncada, procurar-se-á realizar, em Março de 1981, um número da «Técnica» em homenagem ao falecido Professor. 
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* NOTAS CIENTÍFICAS é uma secção recentemente criada na «Técnica», destinada à publicação de resultados 
originais de investigação teórica ou experimental. Os trabalhos remetidos a esta secção deverão ser redigidos prefe- 
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abstract, figuras e referências bibliográficas. 

A Direcção da «Técnica» 


CURSOS NO AMBITO DA GEOTECNIA 


No seguimento dos cursos que anualmente vem organizando desde o ano lectivo de 1975/76, vai a Univer- 
sidade Nova de Lisboa, iniciar em 7 de Outubro de 1980, 2 cursos de pós-graduação, um em Geologia de Enge- 
nharia, especialmente destinados a licenciados em Geologia e em Engenharia de Minas e o outro em Mecânica 
dos Solos, especialmente destinado a licenciados em Engenharia Civil. 

Os interessados poderão obter mais informações nos Serviços Académicos da Faculdade de Ciências e Tecno- 
logia da UNL/Seminário dos Olivais (Moscavide), Lisboa. 
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NOTA CIENTÍFICA N.º 3 


Media where the coincidence of the 

directions of phase and group velo- 

city remain under a Lorentz Trans- 
formation 


M. ABREU FARO 

MARIA JOÃO MARQUES MARTINS 
Universidade Técnica de Lisboa (1.5. T.) 
Centro de Electrodinâmica (Il. N. 1. C.) 


SUMMARY 


E 


c2 
—k.k = C, where C, is a relativistic invariant, the direc- 
tions of the phase and group velocities are coincident 
and remain coincident under a Lorentz Transformation, 
We prove in this paper the converse: Imposing the 
coincidence as a relativistic propriety, the dispersion 
— «g 
equation is necessarily of the form rh k= Cro 


H the dispersion relation is of the for 


SUMÁRIO 


Se a equação de dispersão de um meio for da forma 
—9 


= —k.k = C, em que C, é um invariante relativista, 
2 


as direcções das velocidades de fase e de grupo são 
coincidentes e permanecem coincidentes sob uma Trans- 
formação de Lorentz. 

Provámos o reciproco: Impondo a coincidência consi- 


derada como uma propriedade relativista a equação de 
= 
—k.k=C,. 


dispersão é necessariamente da forma : 
c 


For a homogeneous isotropic medium, in the rest 
frame, the directions of the phase and group velocity are 
coincident. 

Under a Lorentz transformation these directions are 
in general no longer coincident. This is well known, 

Once in an isotropic collisionless plasma the disper- 
sion equation is an invariant, the direction of the phase 
and group velocity remain coincident under the Lorentz 
transformation. 

We may then put the question: 

— What is the dispersion relation of media where we 
could eventually observe the same phenomena? 

We introduce one single condition: 

— |f the directions of the phase and group velocity 
coincide in any frame, then they must remain coincident 
in any other inertial frame, 


Manuscrito recebido para publicação em 18/6/80. 
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Let us consider the expressions of the phase and 
group velocity in the rest frame S(x'), moving with 
velocity v relative to S (x): 


EE Ed e Sit fem = 
Vp=—— k=-—— (ke, +k,e,+k,e,) (1.8) 
k.k k.k , ca fa 
“Rs 
= do  dIw dt) à tw 
V; = = e, + e. + e, (1.b) 
dk Dk, dk, dk, 


IF V, and Va are paralell, then we have 


(2) 


Transforming these quantities to frame S(x!') and 
imposing the coincidence of the directions of phase and 
group velocity, we have 


It is well known that the group velocity is transformed 
in the same way as the velocity of a material particle, so 


dt 
E ——— —w 
Ú tu Ee 0 ks 
iii A 
dk, | V Eidos 
c? dk 
É ty 
dm 1 à k 
mi ES : à tu) (8) 
dk á à ass 
7 c? Ok 
É ty 
dw 1 dk, 
E o a d 
ok, Y E A ci 
c? Ok, 


Once [k, tww/c] is a four-vector, we also have 


= t|——o | k=k k=k (5) 
Li 


Inserting (4) and (5) in (2) we obtain, once we always 
v dm 
have (1 pe ori; o ) = 0 
ct Ri 


dk 
REM 
dk 1 do 1 20 
la RR A ASS ES (6) 
v k, ok, k, Ok, 
K— — 
E: 
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We immediatly verify that the second equality of (3) 
is a consequence of the second equality of (2). 
From the first equality of (3) we obtain 


ú | da 

cmi (o) E (7) 

dk, é k, dk, 

which leads to 
à 2 
dke kw 
According to (3) 

d 2 
My Dc (8.b) 
dk, tu 
A ia ç] 
ds (8.0) 
dk, k ty 


If we take the scalar product of the phase and group 
velocity we have 


dm tu do uy) dt 
ViMp= Dk + k, 42 D+k (9a 
PÉ dk kk CU dk kk dk, kk 


x 


Taking 8.8) b) c) into account we have 
Me =€ (9.b) 


Since we have considered that the directions of 
phase and group velocity are the same, we conclude 


VV. =? (10) 


According to the Principle of Relativity it is obvious that 
from (3) we would also conclude 


dw 1 2 
ssa — (11.,a) 
dk. k. ty 
do 1 2 | 
— — = É (11.6) 
à kk, ty 
RS 2 
dada (11.9) 
dk, k, 1) 
hence 
V, V = c? 


k? + k?) = C, = const. (12) 


This is the dispersion equation of isotropic media 
which satisfy the properties we have postulated, 

This relation is observed, as we have previously 
stated, in an isotropic collisionless plasma where 
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EE 4 
ek k= Es e (13) 


and U, = VA. e? /E,m, is the plasma frequency. 

The first member of (13) is a relativistic invariant and 
therefore 9, must be a relativistic invariant as well. 

As a matter of fact, from the transformation of the 
electronic density N. and the mass m, we can easily 
conclude that to, is an invariant, 

Obviously, equation (12) is also verified when the 
medium considered is the vacuum, where C, = O and (12) 
is then given by 


A -k.k=0 
c? 


Equation (13) is also verified by the formalism intro- 
duced by De Broglie in Wave Mechanics, where 


Moment > p = de 
2% 

h — 

Energy 2 E=—w 
2% 


a 
D =| is a four-vector and so 
E 


RO asim 7 

Ckic Cm 

c2 h? Õ [E] 
CONCLUSION 


If the dispersion equation of a homogeneous isotropic 
medium is 


: — k.k = relativistic Invariant 
E 


then the directions of the phase and group velocity are 
coincident under a Lorentz transformation, 

The purpose of this paper is to demonstrate the 
converse: 


If these directions remain coincident, we necessarily 
have for the dispersion equation 


where C, is a relativistic invariant, and we always have 


ViVe= VV, = 0? 


k is interesting, we believe, to applivy these results 
to a waveguide but this can be developed later on. 
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Determination rapide du courant 
absorbé par un moteur asynchrone 
alimenté par onduleur de tension 


F. LABRIQUE 
J. M. THIAY 
Université Catholique de Louvain 


ABSTRACT 


In this short paper, it is shown that the current drawn 
by an asynchronous motor from a voltage inverter power 
supply can be easily and quickly predicted without resort- 
ing to a harmonic analysis of the inverter output voltage. 


RESUMO 


A presente nota mostra que a corrente absorvida 
por um motor assincrono alimentado por um ondulador 
de tensão pode ser determinada de maneira simples e 
rápida, sem recorrer à análise harmónica da tensão de 
saida do ondulador. 


1 — INTRODUCTION 


Le calcul du courant absorvé en régime permanent 
par un moteur asynchrone triphasé alimenté par un ondu- 
leur de tension s'effectue habituellement en décomposant 
en série de Fourier la tension de sortie de I'onduleur et 
er appliquant chaque composante au circuit équivalent 
de la machine correspondant à l|'harmonique considé- 
rée [1], [2]. En suivant cette procédure, on introduit 
généralement une simplification qui consiste à réduire, 
pour toutes les harmoniques, le schéma équivalent de 
la machine à sa self de fuites. 

Toutefois, lorsque cette simplification est utilisée, il 
est inutile de recourir à un calcul, composante par com- 
posante, pour déterminer le courant harmonique. En effet, 
dans ce cas, on peut calculer directement le courant 
harmonique total à partir de la seule connaissance de 
la tension de sortie de l'onduleur et de sa composante 
fondamentale. On aboutit ainsi à une méthode de calcul 
plus simple et rapide, qui est développée dans la pré- 
sente note, en considérant à titre d'exemple le cas d'une 
alimentation par un onduleur à ondes rectangulaires dont 
la tension côté continu a une amplitude proportionnelle 
à la fréquence de sortie de maniêre à maintenir constant 
le flux dans le moteur (fig. 1) [3]. 


Manuscrito recebido para publicação em 1/7/80. 
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Figure 1b 


Des essais effectuês sur une machine de laboratoire 
montrent la bonne concordance existant entre les pré- 
dictions théoriques et les résultats expérimentaux. 


2 — CALCUL DE LA COMPOSANTE HARMONIQUE TO- 
TALE DU COURANT 


La tension de phase dalimentation, fournie par 
londuleur, est une fonction en échelons dont la forme 
est donnée à la figure 1b pour I'onduleur à ondes rec- 
tangulaires pris comme exemple, 

On prend comme schéma équivalent d'une phase de 
la machine, le circuit représenté à la figure 2 [4]. 


Les composantes harmoniques présentes dans la 
tension d'alimentation sont de rang k=6jÉ1, (je (1,2,3...)) 
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et le glissement du rotor par rapport aux harmoniques 
est donné par l'expression (1) 


(6) = 1) ww, 6) 1 (1) 
"em Ena  g=1. 


Dans cette expression, tw représente la pulsation de 
la fondamentale et tu, la vitesse angulaire de rotation du 
rotor qui, en fonctionnement normal, est sensiblement 
égale à tw. 

Comme le glissement du rotor par rapport aux har- 
moniques est pratiquement égal à 1, on peut, sauf aux 
três basses fréquences d'alimentation ou les chutes de 
tension sur les résistances AR, et A; ne peuvent plus 
être négligées, reduire le schéma équivalent, pour les 
harmoniques, à l'inductance totale de fuites, Lg; (figure 3). 


Fig. 3 


Des lors, la composante harmonique totale du cou- 
rant d'alimentation s'obtient de maniêre immédiate en 
faisant la différence entre la réponse en régime de 
inductance totale de fuites Lf à la tension de phase de 
onduleur, soit |, et sa réponse en régime à la com- 
posante fondamentale de cette tension, soit |'. La compo- 
sante fondamentale de la tension est donnée par |'ex- 
pression (2) 


2E 
Ve= sin mt (2) 
Tt 


Comme la réponse en courant d'une inductance à 
une onde de tension s obtient par simple intégration, les 
résultats de la figure 4 sont évidents en n'appellent aucun 
commentaire. 

La position et l'amplitude des pointes de courant 
harmonique sont reprises au tableau 1. 


COURANT HARMONIQUE TOTAL 


Position Valeur absolue des extrêmés 


ZE 


Tableau 1 
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Figure 4 


3 — CALCUL DU COURANT TOTAL ABSORBE PAR LE 
MOTEUR 


Le courant de phase total absorbé par la machine 
sobtient en ajoutant au courant harmonique total, le 
courant fondamental absorbé soit If. Ce courant se cal- 
cule à partir du schéma équivalent de la machine pour 
la composante fondamentale (figure 5). 


We L 
most 


ERES 


Le courant fondamental peut se décomposer en ses 
composantes active et réactive (l, et Ip) 


V. Rar 
V E Vil! 
wgL, ( Rr y+ 212 


ln = Ipcosg = 


Ip = lesing = 
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Le courant total absorbé par le moteur est repré- 
senté à la figure 6. 


Figure 6 


Les pics de courant principaux sont localisés en 
wt = km et wt=27/3+ kr. Si on admet, en premiêére 
approximation, que la composante réactive du courant reste, 
pour tous les régimes de fonctionnement, pratiquement 
égale au courant magnétisant ut les pointes de courant, 


exprimées en valeur de crête, peuvent s'écrire en fonction 
des grandeurs nominales de la machine Iy et cos by et 
du rapport w Ltu/w,L,, comme valant: 


à vide 


toL 


[Gm |= | 0090. =— É V2 x sin Gn 
toe Le 


ho 


top 


[1 + km) - 0,5 + 0,048 E) V2Ix sin by 
| to Le 


( 2% 
3 
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à la charge nominale 


| W, a = 
| | (km) - 1 + 0.0966——— | V2Iysin ón 
Lg Le 
| (— + km) | = 
weL 6 
=| 0,5 + 0,048 n/2 Isin + ly COS É 
| PÇA V2Iysin ou E. N 


4 — VERIFICATION EXPERIMENTALE 


Les essais ont été réalisés sur un moteur asyn- 
chrone à rotor bobiné d'une puissance de 7,5kW dont 
les paramétres à la fréquence d'alimentation utilisés sont: 


wW L = 21.420 
Eu 
A) = 0910 


L'onduleur fonctionnait avec une tension E de 114,8V 
correspondant à une fréquence d'utilisation de 23,8 Hz. 

Les ondes de courant relatives à 2 régimes diffé- 
rents de la machine ont été relevées. Elles correspon- 
dent à des glissements de 0,85% et 12,4%. 

Les figures 7 et 8 présentent les courbes relevées 
expérimentalement ainsi que les valeurs déterminées par 
calcul et montrent la três bonne concordance existant 
entre la théorie et |'expérience. 


5 — CONCLUSION 


On a présenté une méthode de caleul simple et 
rapide du courant absorbé par un moteur asynchrone 
alimenté par un onduleur de tension à ondes rectangu- 
laires. 

Cette méthode, dont la précision a été vérifiée expé- 
rimentalement, peut s étendre sans difficulté aux cas 
d'une alimentation par un onduleur à modulation de lar- 
geur d'impulsion pour lequel le recours à la série de 
Fourier serait encore plus compliquée que dans le cas 
de |'onduleur à ondes rectangulaires qui a été considéré 
dans la présente note. 
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Resumo dos Artigos publicados na «Técnica» n.º 459 
Ano LV — Junho 1980 


FRANCISCO CARVALHAL, ISABEL LOURTIE e 
PEDRO LOURTIE — Projecto apoiado em computa- 
dor de sistemas de controlo, «Técnica», Lisboa, LIV 
(459) Junho 1980, p. 355-368. 


O artigo descreve com algum detalhe o bloco opera- 
cional de um conjunto articulado de programas para o 
projecto apoiado em computador de sistemas de con- 
trolo. Neste bloco encontram-se implementados os mé- 
todos do lugar geométrico das raizes, Nyquist, Bode 
e da simulação digital da resposta no tempo. Apresen- 
ta-se os algoritmos utilizados bem como alguns exem- 
plos de aplicação. 


MARIA CELESTE ANTUNES e SILVESTRE DIAS 
ANTUNES — Tecnologia da metanogénese (revisão 
descritiva), «Técnica», Lisboa, LIV (459) Junho 
1980, p. 369-382. 


Neste artigo caracterizam-se as variáveis actual- 
mente consideradas como fundamentais no processo me- 
tanogenético, descrevem-se as diferentes arquitecturas 
mais correntes nos digestores rurais, industriais e labo- 
ratoriais, sem que se faça qualquer análise quantitativa 
das relações matéria-prima — produtos. No entanto pro- 
cura-se tanto quanto possível apresentar uma bibliogra- 
Ra exaustiva do estado actual da tecnologia da metano- 
génese, 


SILVESTRE DIAS ANTUNES e MICHEL J. INS- 
TALLE — Utilização do plano de fases no controlo 
da digestão anaeróbia (aplicação à metanogénese), 
«Técnica», Lisboa, LIV (459), Junho 1980, p. 383- 
-392. 


Estuda-se neste artigo a utilização da análise no 
plano de fases no controlo da actividade biológica de 
digestores homogéneos, sem que se faça a lavagem da 
biomassa presente, ou se entre em sestado de choques. 


F. CALDEIRA SARAIVA — Financiamento por acções 
de preço flutuante, «Técnica», Lisboa, LIV (459), 
Junho 1980, p. 393-402. 


O presente artigo é uma aplicação de técnicas de 
controlo à determinação da política óptima de financia- 
mento para firmas cujas accões tem valores de mercado 
decrescentes com a taxa de endividamento. Por aplica- 
ção do Princípio de Pontryagin em tempo continuo as 
políticas dinâmicas óptimas são obtidas como funções 
o rarores dos parâmetros, sendo os resultados finais 

iscutidos. 


| 


J. R. CALDAS PINTO —— Simulador industrial de 
controlo de nivel, «Técnicas, Lisboa, LIV (459), 
Junho 1980, p. 403-408. 


Para fins didáticos construiu-se em 1976/77, no Cen- 
tro de Automática da Universidade Técnica de Lisboa 
(C. A. U,. T. L.), um sistema para estudo do controlo 
de nível, utilizando sempre que possível componentes 
industriais. Na nota técnica anexa descreve-se suma- 
riamente este sistema. 


M. ABREU FARO e MARIA JOÃO MARQUES MAR- 
TINS — Meios em que as direcções das velocida- 
dades de fase e de grupo coincidem e se mantêm 
coincidentes sob uma transformação de Lorentz, 
«Técnica», Lisboa, LIV (459), Junho 1980, p. 409- 
-410. 


Be a equação de dispersão de um meio for da forma 


= di k.k=C, em que C| é um invariante relati- 
vista, as direcções das velocidades de fase e de grupo 
são coincidentes e permanecem coincidentes sob uma 
Transformação de Lorentz. 

Provámos o recíproco: Impondo a coincidência con- 
gsiderada como uma propriedade relativista a equação de 


dispersão é necessariamente da forma = K. k=C,. 


to 
o? 


F. LABRIQUE e J. M. THIRY — Determinação rápi- 
da da corrente absorvida por um motor assincrono 
alimentado por um ondulador de tensão, «Técnicas, 
Lisboa, LIV (459), Junho 1980, p. 411-412. 


A presente nota mostra que a corrente absorvida 
por um motor assíncrono alimentado por um ondulador 
de tensão pode ser determinada de maneira simples e 
rápida, sem recorrer à análise harmónica da tensão de 
saida do ondulador. 
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J. R. CALDAS PINTO — Industrial simulador for 
level control, «Técnica», Lisboa, LIV (459), Junho 
1980, p. 403-408. 


For didactical purposes, a level control system has 
been assembled at CA U. T. L., in 1976/77. Industrial 
components were preferably used. This technical paper 
presents a brief description of the system. 


M. ABREU FARO e MARIA JOÃO MARQUES MAR- 
TINS — Media where the coincidence of the direc- 
tions of phase and group velocity remain under a 
Lorentz Transformation, «Técnicas, Lisboa, LIV 
(459), Junho 1980, p. 409-410. 


If the dispersion relation is of the form 


— K.k = C where C is a relativistic invariant, the 
directions of the phase and group velocities are coinci- 
dent and remain coincident under a Lorentz Transior- 
mation. 

We prove in this paper the converse: Imposing the 
coincidence as a relativistic propriety, ua dispersion 


equation is necessarily of the form = es = Co. 


F. LABRIQUE e J. M. THIRY — Quick calculation 
of the current in an induction motor driven by a vol- 
tage inverter, «Técnica», Lisboa, LIV (459), Junho 
1980, p. 411-412. 


In this short paper, it is shown that the current 
drawn by an asynchronous motor from a voltage inver- 
ter power supply can be easily and quickly predicted 
without resorting to a harmonic analysis of the inver- 
ter output voltage. 
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FRANCISCO CARVALHAL, ISABEL LOURTIE e 
PEDRO LOURTIE — Computer-aided control system 
design, «Técnica», Lisboa, LIV (459) Junho 1980, 
p. 355-368, 


The paper describes with some detail the operatio- 
nal block of a control systems computer-aided design 
package. The root locus, Nyquist. Bode and time res- 
ponse methods are implemented in this block. The algo- 
Fra are presented as well as some application exam- 
ples. 


MARIA CELESTE ANTUNES e SILVESTRE DIAS 
ANTUNES — Anaerobic digestion technology (re- 
view), «Técnica», Lisboa, LIV (459), Junho 1980, 
p. 369-382. 


In this artície, variables nowadays considered as 
fundamental in the anaerobic digestion process are cha- 
racterized, and the different architectures of digesters 
used in rural, industrial and laboratory applications 
are described. Although, if quantitative analysis of the 
transformation: prime matter — products is not consi- 
dered, an exhaustive bibliography on the present state 
of the anaerobic digestion technology has been tried, 


SILVESTRE DIAS ANTUNES e MICHEL J, INS- 
TALLE — Anaerobic digestion control using the 
phase-plane analysis [with application to methano- 
genesis), «Técnicas», Lisboa, LIV (459), Junho 1980, 
p. 383-392. 


Phase-plane analysis of homogeneous digester sys- 
tems is used to give information on how to control the 
biological activity without washing out the reactor or 
reaching a eshock states. 


F. CALDEIRA SARAIVA — Corporation financing 
under falling share prices, «Técnica», Lisboa, LIV 
(459), Junho 1980, p. 393-402. 


The present paper is an application of control tech- 
niques to the determination of the optimal financial 
policy, for a firm whose shares" market value falls 
with increasing leverage. By the application of Pon- 
tryagin's principle in continuous time the optimal dyna- 
mic polícies are derived as funtions of the values of 
the parameters and the results are discussed, 
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PEDIDOS À [ERNIES 


If you could talk to Thomas Edison, 
he'd tell you what it was like to turn the 
lights on in 1879. You could tell him about 
some bright ideas of the 20th century... 
particularly, a technological phenomenon 
that can handle everything from solar heat 
control to lighting your home via voice 
command. The Apple personal computer. 


Expand yourown . 
inventíveness with the 
always-expandable Apple. 

Take a look inside your local computer 
store. There's a range of Apple systems 
for you...whether you want expansion 
capabilities of four or eight accessory 
slots...or memory expandable to 64K 
bytes or 128K bytes. With this kind of flexi- 
bility, the possibilities for creating your 
own computer system are endless. 

Want to add an A to D conversion 
board? Apple makes it happen. Want to 
plug into time sharing, news and elec- 


TO FACAMENT 


With Apple, Edison could've written a 
program to determine why some filaments 
burned longer than others. 


tronic mail services? Apple does it all. 
Because Apple is the most popular per- 
sonal computer with the least complicated 
interface, over 100 companies supply 
peripherals for the Apple family...includ- 
ing an IEEE 488 bus for instant control. 


Disk dríves, a tool kit 

and creativity ín color. 
Apple was one of the first to use 

disk drives for increased performance 

and application versatility. Today, our 54º: 

disk drive offers high density (143K bytesh 


high speed and low cost. No 
wonder this drive is the most 
popular on the market. 

But now Apple goes one 
better with the DOS Tool Kit. 
A series of utility programs, 
it gives you the freedom to 
easily design 280hx192v 
graphic displays in a palette 
of living color... depending 
on your choice of Apple 
system. 

Edison was first with the 
movie camera and projector. 
Now, with Apple's DOS Tool 
Kit, you can be first to work 
wonders with colorful 
creative animation. 


Imagine the 
broadest line 
programs ever. 

Apple's broad line of ofanimation. 
peripherals is equalled only — “ 
by the most extensive line of software 
you'll find in the personal computing 
world. Since more than 170 companies 
offer software for the Apple family, you 
can have one of the most impressive 
program libraries ever. 

When you write your own programs, 
your Apple speaks creatively in BASIC, 


ésappie computer 


Edison had the first movie 
camera... and Apple has 
the DOS Tool Kit that takes 
you into the colorful world 


6502 assembly language. 
Use these languages to score 
a sonata. Apple will play 
back your musical master- 
piece on its built-in speaker. 

Edison listened to his 
voice on a revolutionary 
phonograph in the 1800s... 
now you can listen to the 
sounds of today with Apple's 
inventive family of personal 
computers. 


| 
| 
| Pascal, FORTRAN, PILOT and 


Where to find even 

more illuminating 

Apple experiences. 
There's always something 

new being invented at Apple to 

set your imagination soaring. 

And there's always an expert 

to tell you all about it in 

detail. Your Apple dealer. If you 

already own an Apple, there's 

a whole future ahead to 


challenge man, mind and machine. 


Ifyou're considering a personal com- 
puter, stop by the computer store and 
compare. Apple's reliability, proven perfor- 
mance and recognized technological leader- 
ship will help you see the light. Don't 
let history pass you by. 

Visit our sole distributor 

in Portugal, Sorubal S.A.R.L,, 
Rua Gen. Pimenta de Castro 
lô 8.º, 1700 Lisboa, 

or call 8965 55 (Lisbon) 
telex 12775 SORBAL P, 


